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1.  INTRODUCTION 

The  problem  of  diffraction  of  plane  waves  through  a slit  in  a 
perfect  electric  conductor  of  finite  thickness  has  been  studied  by 
several  investigators  [1-4].  The  most  extensive  investigation  was  that 
of  Lehman  [1],  who  used  the  analytic  properties  of  finite  Fourier  trans- 
forms. The  solution  of  Kashyap  and  Hamid  [2]  used  a Wiener-Hopf  and 
generalized  matrix  technique.  Both  of  these  solutions  were  done  for  the 
TM  case  (incident  electric  field  parallel  to  slit  axis).  The  solutions 
of  Hongo  [3]  and  of  Neerhoff  and  Mur  [4],  were  obtained  bv  a numerical 
solution  of  coupled  integral  equations  and  were  done  for  the  TE  case. 

In  this  report,  we  use  the  method  of  moments  to  solve  coupled  integral 
equations  similar  in  form  to  those  derived  in  [4]. 

This  report  utilizes  the  generalized  network  formulation  of  coupling 
through  apertures  developed  in  [5]  and  [b]  and  extends  these  results  to 
three  regions  coupled  by  two  apertures.  To  accomplish  this  the  equiva- 
lence principle  is  used  to  replace  both  faces  of  the  slit  by  perfect 
conductors,  each  of  which  carry  magnetic  current  sheets  on  both  sides. 

The  original  problem  is  now  broken  up  into  three  regions  which  are 
coupled  by  the  postulated  magnetic  current  sheets.  The  two  half  space 
regions  are  loss-free  with  arbitrary  u and  e and  the  medium  in  the  slit  is 
assumed  lossy  with  arbitrary  complex  u and  r. 

Continuity  of  the  tangential  magnetic  field  is  used  to  derive  two 
coupled  operator  equations  involving  the  equivalent  magnetic  currents  as 
unknowns.  These  equations  are  pit  into  matrix  form  using  the  method  of 
moments,  and  solved  by  using  standard  matrix  methods. 


I 


The  aperture  coupling  between  the  three  regions  is  characterized 
hv  a combination  of  "admittance  matrices"  computed  separately  for  each 
region.  This  gives  rise  to  a network  interpretation  of  the  problem 
which  treats  the  unknown  magnetic  currents  as  port  voltages  and  the  ex- 
citation as  port  currents. 


1 Li.  TKPBl-KM  kormvi.ation 

The  original  problem  configuration  is  shown  in  Fig.  1.  It  con- 
sists of  a perfect  electric  conductor  of  thickness  d separating  two 
regions  a and  c which  may  have  different  electrical  properties.  Coupling 
between  the  two  regions  occurs  through  a slit  of  width  w filled  with  «n 
arbitrarily  lossy  medium.  The  conductor  is  infinite  in  the  z and  y 
directions.  The  problem  consists  of  three  regions  separated  by  two 
boundaries  (the  slit  faces)*  Using  the  equivalence  principle,  the  three 
regions  can  be  separated  by  covering  the  slits  with  perfect  electric  con- 
ductors and  magnetic  currents,  as  described  in  [5]. 

The  regions  are  defined  by 


region  a 
region  b 
re  g i on  c 

and  the  boundaries  as: 

T* 

' 1 

‘ 2 

which  are  the  two  boundaries  of  s 
principle,  !' ^ and  I'  are  covered 


x < 0 , ally 

0 < x < d , 0 < y < w 

x '*  d , all  y 


X 

N 

O 

0 < 

y 

< w 

X - d. 

0 < 

V 

< w 

'parat ion.  To  utilize  the  equivalence 
'v  perfect  electric  conductors,  and  on 


each  side  of  these  conductors  a magnetic  current  sheet  is  placed  which 


1h  determined  by 


Nj  - li  * Kj  on  l'j  (la) 

- K.,  * ri  on  I’2  (lb) 

where  F^  and  are  the  total  electric  fields  in  the  aperture  at  I'j  and 
r2  respectively.  This  breaks  the  original  problem  up  into  three  parts 
as  shown  in  Fig.  2.  The  magnetic  current  sheets  on  each  side  of  F.  and 
T.,  are  equal  and  opposite  In  sign  for  each  slit  because  in  the  actual 
problem  the  electric  field  must  be  continuous  across  them. 

The  tangential  component  of  magnetic  field  In  region  a at  boundary 
I'j,  H't*j,  is  due  to  the  Incident  field,  H*,  and  also  to  the  magnetic  cur- 
rents Mj  radiating  in  the  presence  of  a complete  perfectly  conducting 
plane  at  x - 0.  This  is  written  in  operator  notation  as 


ii“,  - »“,(«,>  + iitl 


(2) 


— mm  — | 

where  ll(j(Mj)  iltu*  Hj.  Rive  the  field  of  sources  radiating  In  the  presence 
of  the  complete  conductor,  and  subscript  1 is  used  to  denote  that  the 
field  point  is  on  I'j. 

In  region  b,  where  there  are  no  impressed  sources,  the  tangential 
component  of  magnetic  field  on  I'j  is  due  to  the  equivalent  currents  -Mj 
and  -M2  or 

"u  - ♦"V-v 


- H'tVV  - 


(3) 


where  the  minus  signs  factor  out  because  the  field  operatoi s are  linear. 
Also  in  region  b the  tangential  magnetic  field  at  l\,  is  written  as 


Fig.  2.  Equivalences  for  regions  a,  b,  and  c. 


nt2  ' fit2<-V  + nt2<-«2> 

• - H^Fy  - fi^2(M2)  (4) 

where  subscript  2 indicates  that  the  field  point  is  on  ?2 . 

In  region  c,  where  there  are  also  no  Impressed  sources,  the  tangential 
component  of  magnetic  field  at  I’.,  is  due  only  to  currents  M0  or 

fit2  ■ ^2(S2>  <5) 

— c — 

where  H gives  the  field  of  sources  radiating  in  the  presence  of  a 

complete  conductor. 

The  true  solution  to  the  original  problem  is  obtained  when  the 
tangential  components  of  magnetic  field  are  continuous  on  the  two  boundaries 
Tx  and  I',,.  Thus  setting  (2)  equal  to  (3)  and  (4)  equal  to  (5)  we  have 

nt\(V  + ”tl(V  + "tV*^  = ~ (6a) 


Ht2(V  + "t2(M2)  + Ht2(M2)  * ° 


These  two  coupled  operator  equations  must  be  solved  simultaneously  and  we 
apply  the  moment  method. 

The  first  step  is  to  assume  that  the  two  magnetic  current  sheets  may  be 
expanded  by  a linear  combination  of  basis  functions  defined  on  and  1.,  as 


M = y V,  M,  on  T, 

1 In  In  1 


S2  ’ l V2„52„  r2 


Here  V and  V.  are  unknown  complex  scalars  and  M.  and  M_  are  vector 
In  2n  r In  2n 

basis  functions  on  T ^ and  I",,  respectively. 
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~ - - - — ■— a — • - ■ m — — — - — - — ^ — 


Now  Hub at  it u ting  Kq».  (7a, b)  Into  Kqa.  (ba,b)  we  obtain 


"l  “l  N2 

t V'ci'V  + \ VlVV  * l,  V"-.'",.'  - - 

n-1  n-1  n-1 


(8a) 


Nl 

V 

L 

n-1 


i , * » * 

l V,  H®  (M.  ) + y V.,  H®  (M_  ) + V V,  fi*\  (M  ) - 0 
L , in  t2  In  *•,  2n  t2  7n  L.  2n  t/  2n 
n-l  n- 1 


(8b) 


An  appropriate  symmetric  product  Is  defined  as 


<X,  Y>  - 


X • Y dy 


ri*  r2 


(<») 


where  the  variable  of  Integration  Is  either  on  F^  or  • Also 

needed  are  a set  of  testing  functions  {W^,  m - 1,2,...,  N^)  on  Tj  and 

{W.,m,  m • 1,2,...,  N.,)  on  F,.  Next  take  the  symmetric  product  of 

(8a)  with  W.  and  (8b)  with  W,  to  obtain 
lm 

N . 

I V,  ^Wi  )>  + <W.,H°  (M  M 

L In  lm  tl  In  lm  tl  in 


n-1 


+ X V-  <W1  .H?1  (M,  )>«-<W  ,\'}> 

L . ^n  lm  tl  2n  lm  t 

n-1 


(10a) 


for  m-  1,2, ...  ,N^,  and 
N1  N2 

X vln  'W2m-  Ht2(Mln)>  + * v*«  l <**«>>  + <W,m.  <2(M,n)>] 


n-1 


n-1 


2n 1 2m  1 2 2n 


2m*  t2  in 

- 0 (10b) 


for  m - 1,2,...,  N,.  These  equations  raav  be  rewritten  In  matrix 
form  as 
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[Y21]V1  + [Y22  + Yh8C)V2  - 0 

where  the  various  component  matrices  are  explicitly  identified  by 


(lib) 

| I 


[Yh8a]  - - 

(12a) 

[Yh8C]  = - 

<2(M2„)>'n2.N2 

02b) 

[YU]  - - 

l<Wln' 

(12c) 

[Y12 1 - - 

(12d) 

[Y21 1 - - 

Ht2(Mln)>]N^  Nj 

( 12e) 

1 1 

1 

[Y22]  = - 

Ht2(M2n)r>1N2x  N, 

( 1 2 f ) 

i<ww 


( 12  r) 


Kquations  (11a)  and  (lib)  compose  a (N^  + ^)  * (N^  + N2)  system  which 
suggests  the  network  representation  shown  in  Fig.  3 where 

[Y11]  [Y12]' 

[Yb]  - (1 

[Y21]  [Y22] 

tl 8 <1  V)  j^gc . 

The  matrices  [Y  ],  [Y  ],  and  [Y  ] are  the  network  representations 


of  regions  a,  b,  and  c respectively.  The  explicit  computation  of 
these  quantities  which,  to  carry  the  network  analogy  further  we  call 
admittance  matrices,  depends  only  upon  their  respective  regions. 
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HI.  PROBLEM  SPECIALIZATION 

The  development  up  till  now  has  been  quite  general.  In  this 
section  we  define  the  basis  functions  used  in  the  expansion  and  test- 
ing procedure.  Also  the  electrical  properties  of  the  three  regions  are 
described. 

A z-directed  magnetic  current  excitation  produces  a field  TE  to 
z,  which  requires  only  z-directed  equivalent  magnetic  currents  over 
V1  and  r2<  The  basis  functions  on  are  chosen  as 


1 u^  , (n-1)  Ay  £ y £ nAy 


In 


(14) 


0 ’ 


elsewhere 


for  n-l,2,...,N  and  Ay  » — , y being  defined  on  T . M is  defined 

‘N*  1 2n 

exactly  the  same  but  with  y defined  on  T^.  Here  we  have  taken  - N 

so  that  the  number  of  expansion  functions  on  each  slit  face  is  the  same. 
The  testing  functions  are  chosen  identical  to  the  expansion  functions  so 
that  and  - Mj  . Thus  from  Eqs.  (12c  through  f)  we  have 


[YU]  - [Y11  ] 


[Y21]  - [Y12 ] 


(15) 


[Y22]  - [Y22] 


where  the  tilda  (')  denotes  the  usual  matrix  transpose.  Also  since  the 
basic  functions  on  r2  are  the  same  as  those  on  Tj,  we  have 


11  22 
lY11]  - [Y“] 


[Y21]  - [Y12] 


(16) 


as  will  be  seen  more  explicitly  in  Section  TV-3. 
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The  electrical  properties  of  regions  a and  c are  specified  by  the 

real  numbers  e /e  , p /p  , e /e  and  p /p  where  c and  p denote 
aoaoco  co  o o 

permittivity  and  permeability  of  free  space.  The  properties  of  region 
b are  given  by  the  complex  numbers  c^/eo  ®nd  p^/pq  to  describe  the  losses 
present . 


IV.  COMPUTATION  OF  MATRIX  ELEMENTS 

lis  d 

1.  Admittance  matrix  for  Region  a.  [Y  ].  These  matrix  elements 
are  found  by  computing  Eq.  (12a).  Since  H^(M.^)  gives  the  field  of  cur- 
rent sheet  M,  radiating  into  the  half  space  a in  the  presence  of  a com- 
ln 

plete  conductor,  we  may  write 


y-y' | )dy ' 


(17) 


Here  Ay  is  the  region  on  T,  where  M,  i*  0,  k “w  /p  e , r\  ■ rp It  , 

n 1 In  a aaaaa 

(2) 

u)  ■ radian  frequency  of  line  source,  and  is  the  Hankel  function  of 

second  kind,  order  zero.  Substititing  Eq.  (14)  for  M.  ani'  using  the 

fact  that  W,  - M,  for  m-l,2,...,N,  we  have  for  the  mnth  element  of 
ltn  lm 

Eq.  (12a): 


yhsa 

mn 


IT  \ l H<2>  (kJy-y'IWy  dy’ 
“ Aym  Ayn 


(18) 


a 

From  this  equation  it  is  evident  that  [Y  ] is  a symmetric  Toeplitz 
matrix,  i.e.,  the  elements  are  functions  onlv  of  |m-n|  and 
hence  only  one  column  need  be  computed. 
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2.  Admittance  matrix  for  Region  c. fY 


The  elements  for  the  admit- 


tance matrix  of  region  c are  found  by  computing  F.q.  (12b).  The  operator 
**t2^^2n^  8lves  the  of  a magnetic  current  sheet  radiating  into 

region  c in  the  presence  of  a complete  perfect  conductor.  This  is  written 


ir,(K,  ) - - -i-  M,  H'  '(k  y - y’ 
t2  2n  2n  2n  o c 


where  Ay  is  the  region  on  T»  where  M~  i 0,  k “ u>/n  c , and  n - *n  /r  . 

'n  2 2n  ccc  ccc 

Substituting  Eq.  (19)  into  Eq.  (12b)  and  using  the  fact  that  M.,  - W,  we 

•.  n * n 

have  for  the  mnth  element  of  Eq.  (12b): 


- f f Ho2'(kc'y-y,')d*  dy* 


Ay  Ay 
m n 


(Y  L ] is  also  a symmetric  Toeplitz  matrix. 


3.  Admittance  matrix  for  Region  b,  [Y^  ] . The  elements  of  [Y^ ] are 
fcund  by  computing  Eqs.  (12c  through  f) . The  operators  or  H^,  give  the 
fields  of  current  sheets  or  M.,^  radiating  inside  a closed  conducting 
box  as  shown  in  Fig.  4.  Thus  bruaklng  the  computation  of  [Y^l  up  into  two 

parts  we  have 

1121 

Part  l:  Source  M on  1'^,  computation  of  [Y  ] and  [Y  ]. 

- 12  22 

Part  11:  Source  M,  on  !*«,  computation  of  (Y  “ ] and  [Y 

«.  n a. 

For  the  TE  case  under  consideration,  all  the  magnetic  currents  are 
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rectangular 
perfectly 
conducting  box  B 


Part  I 


rectangular 
perfectly 
conducting  box  B 


Part  II 


Fig.  4.  Region  b broken  up  into  two  problems  each  with  an 
equivalent  source. 


■ 

i 


z-dlrected.  Thu*  an  electric  vector  potential  function  F can  be  de- 
fined in  region  b as  [7,  sec.  3-12]: 


F - *u  (21) 

z 

where  * scalar  wave  potential  satisfying  the  two  dimensional  dif- 

ferential equation 

-Mr  + Mjr  k^  ♦ - 0 (22) 

dx*-  3y 


everywhere  in  region  b except  where  the  sources  are. 
gat  ion  constant  in  region  b,  k^»  is  given  as  k^  " wr'u 
are  then  determined  from  F by: 


The  wave  propa- 


.e..  The  fields 
b b 


E 


u T*-  ♦ U T* 
x 3y  y 3x 


(23) 


H - - (24) 

The  tangential  component  of  K obtained  from  Kq.  (23)  evaluated 
on  the  box  B in  Fig.  4 must  be  zero  except  where  there  are  magnetic 
surface  currents.  Hence  for  Part  l,  Fq.  (22)  must  be  solved  in  region  b 
subject  to  the  conditions: 


and 


(25) 


8x 


(26) 
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r 


which  are  obtained  from  the  components  of  Eq.  (23) . The  solution  to  (22) 
satisfying  Eq.  (25)  is  of  the  form 


ill  = T A cos  k (x-d)  cos 

p=0  p xp 


(27) 


where 


k2  - 1?  - <^>2 

xp  D W 


(28) 


The  coefficients  A are  found  by  satisfying  Eq.  (26)  which  is  rewritten 
P 


as 


y A k sin  k d cos  Z-  J~  = - M 

*-  n vr*  vn 


p=0 


p xp  xp 


ZUL 

w “In 


(29) 


Multiplying  both  sides  of  Eq.  (29)  by  cos  and  integrating  from  0 


to  w on  r , we  obtain 


-e 


p w k sin  k d 
r xp  xp 


M cos  dy 
In  w 


-b 


where  e is  Neumann's  number  (e  =1  for  p=0  and  e =2,  p > 0) . Now  H _ 
p o r p tl 


(30) 


and 


r.b 


are  simply  given  by  Eq.  (24). 


Thus  H 1 (M.  ) becomes 
1 1 In 


Hti(V 


- jaicb<»(Mln) 


IWE.  00  e cos  k d f , 

_b  y p..  xp  M COS  dy'  cos  ^ 

w k sin  k d I In  w w 

p=0  xp  xp  ' 


(31) 


with  a similar  result  for  Substituting  into  Eqs.  (12c  and  e)  we 

11  21 

obtain  the  mnth  elements  of  [Y  ] and  [Y  ] : 
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,11 


Jwe  ® c cot  k d 

D " **  ”** 


“ e cot  k 0 r r 

■.-vM  Jv 


cos  cos  £12.  dy  dy’  (32) 

lm  w w ' 7 


r r 
l l 


f21 

mn 


jute.  ® e esc  k d r r 

l^KT^  l 


w *•-  k 

p-0  xp 


Mln  W2mCOS  £?L_  C0S  ^ dy  dy’  (33) 


r r 

2 'i 


For  Part  II  as  shown  in  Fig.  4,  Eq.  (22)  must  be  solved  subject  to 


_3£ 

3y 


3ip 


y»0,  w 


3x 


- 0 


(34) 


x-0 


and 


3ifr 

3x 


M, 


2n 


(35) 


which  are  obtained  from  the  components  of  Eq.  (23).  The  solution  to  (22) 
satisfying  Eq.  (34)  is  given  by 


♦ ■ 7 B cos  k x cos  y- 

p-0  p xp 


(36) 


where  k^  is  again  given  by  Eq.  (28).  The  coefficients  are  found  by 
satisfying  Eq.  (35)  which  is  rewritten  as 


7 - B k sin  k d cos  - M- 

_«  P xp  xp  w 7n 


p-0 


(37) 


Again  multiply  (37)  by  cos  ^^2-  and  integrate  from  0 to  w on  T.  to  get 

w 2 


B - - 
P 


wk 


f 

sin  k d I 


xp  xp 


«2„  coi  dy 


(38) 
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(M.^)  is  again  given  by  Eq.  (24)  j 


\i<V  ' - J“V<V 


Jut. 


w k sin  k d 

p-0  xp  xp 


M„  cos  dy'  cos 

2n  w J w 


(39) 


with  a similar  result  for  Ht2^2n^‘  ^*ua  substituting  into  Eqs.  (12d 

12  22 

and  f ) we  obtain  the  mnth  element  of  [Y  ] and  [Y  ] :> 


f12 

mn 


jue.  00  e esc  k d r 

k Y P XP 


p-0  xp 


W M_  cos  PJ'y. - cos  dy  dy' 
lm  2n  w w 


1 1 2 


(40) 


Jue.  00  e cot  k d 
{22  - _ . > J*  -P  *P 

USX 


w Ln  k 

P-0  xp  r 

2 1 2 


| W2„  «2n 


cos  cos  dy  dy  ’ 

w w 


(41) 


More  is  said  about  the  computation  of  the  admittance  matrix  elements 
of  region  b in  Appendix  A. 


4.  Excitation  Matrix.  The  source  which  is  placed  in  region  a at 

coordinates  (x  , y ) is  a magnetic  current  filament  Ku  radiating  in 
s s z 

tie  presence  of  a complete  conductor  at  x - 0.  Hence  the  tangential 
component  of  incident  magnetic  field,  H*,  is  given  by 


. k K 

hJ  - - u H(2)(k  R ) 

t z 2n  o a s 

a 


(42) 


where  R - /x2  + (y-y  )z.  Substitution  of  this  H*  into  Eq.  (12g)  yields 

S 8 S t 

the  following  formula  for  the  mth  element  of  vector  I*: 
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£ f H<2>OcR)  dy 
In  J o a a 7 

® A .. 


n • 1,2,...,N.  The  variable  of  Integration,  y,  Is  on  1'  and  Ay  la  the 

1 tn 

region  where  W i*  0.  Equatlona  (11a, b)  may  now  be  solved  for  V and  V 


V.  TRANSMISSION  COEFFICIENT 


The  tranamlaalon  coefficient  of  the  silt  Is  defined  as 


trans 


where  ptriinB  1*  the  time  average  power  transmitted  Into  region  c by 

the  ellt  and  P^n{.  Is  the  time  average  Incident  power  Intercepted  by  the 

ellt  from  region  a,  both  for  a unit  length  In  the  t direction.  P 

• trans 

then,  la  Juat  the  real  part  of  the  Poyntlng  vector  flux  through  boundary 


r2  which  la  given  by 


w 

- _*  A 

P.  ■ Re  E0  « H_  • n 
trans  122 


where  * denotea  complex  conjugate.  E^  and  H2  are  the  total  electric  and 
magnetic  fields  at  I^.  Using  the  vector  Identity  [7]: 


E2  x H2  • ri  - H2  • (n  x E2) 


and  Eq.  (lb)  we  obtain 


trara 


Re  »2  • M2dy 


Now  substitute  Eq.  (7b)  Into  the  above  to  get 
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Ke  l V.  { H?  • M„  dy) 


The  magnetic  field  H,  is  due  to  the  magnetic  current  sheet  M,  radiating 


into  region  c or 


t Q 


"j  ' fl;VS2>  - 

n»l 


Substituting  Eq.  (49)  into  (48)  and  using  Kq.  (12b)  and  the  fact  that 

M,  - W,  , we  obtain 
7n  2n 


r - r.  «,[vh"l'Tv* 


This  is  the  usual  formula  for  power  flow  into  the  network  represented 
by  [Yh*C ] of  Fig.  3. 

The  time  average  power  radiated  into  whole  space  per  unit  length 
in  s by  a magnetic  line  source  of  strength  K is  given  by  [7]: 

- *S-  l«l*  <51 

a 

P is  that  portion  of  Pf  which  is  intercepted  by  the  aperture  and  is 


given  by 


P .ip  .Ijl.  |k|2 

inc  2n  f 8irn  11 


where  0 is  defined  in  Fig.  5.  Equation  (44)  is  finally  written  as 


6k  |K 

a1 


Re  V*> 

m I * * l 
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VI  . JVWTK  liAJN^XNJ^MKASj’KI-MKNl  Y1YT0R 

The  power  gain  pattern  1 1\  region  e is  del  ined  as  the  rat  io  ot  the 


radiation  intensity  whieit  would  exist  it  the  transit. itted  time  average 

power  were  radiated  unitormlv  over  halt  spaee  to  1'  in  It; . t**l,  ot 

* 1 1 1 uns 


C(f) 


■nr  n h 

m e m 


t S 


t rans 

is  the  eomponent  ol  the  magnet  le  field  in  region  c in  the  dlreet  ion 


of  a magnet  ie  test  line  current  K n due  to  curia  nt.  sheet  M,  radiating 


in  the  presence  of  a complete  conductor  .it  x - d.  K u is  used  to 


-l 


measure  11  at  position  (r  ,,*■)  hv  reciprceitv.  It  11,.  is  the  field  at 
m ' m } 1 K 


, due  tv'  K u,  radiating  in  the  present  e ot  a complete  conductor,  then 


t t om  i oc i proc  i t v 


K u 

t'l  z 


% 


M,  dv 


Substitution  v't  he; . (?bl  into  (.SSI  results  in 


1 K - J V.,  vM  , II1 
m m * dn  _n  K 
n- 1 


which  is  rewritten  in  matrix  form  as 


11  K - lm  V, 

m V)  1 


071 


1 is  a measurement  vector  defined  as 


0,0  YVi 


08' 


The  elements  ot  1 are  essentially  the  same  as  those  of  T * anil 


are  given  by 


Jl 


i J 


tv*  1 ,2  , . • . ,N . If  r \ (far  field  measurements)  where  \ - ~ then 

me  c k 


VII.  NUMERICAL  EXAMPLES 


All  the  examples  Riven  here  are  done  for  normal  incidence  with 
the  source  far  enough  away  (x  - -100\  ) and  the  strength  adjusted 

S d 

to  simulate  i unit  incident  plane  wave  at  slit  face  Tj  . The  permea- 
bility and  permittivity  of  regions  a,  b,  and  c are  that  of  free  space 
except  where  noted.  To  check  the  computer  program,  magnetic  currents 

were  computed  for  a relatively  thin  slit  (w  • .4\  , d * .001 A ) for 

a a 

different  values  of  permittivity  and  permeability  in  region  c.  These 
results  are  given  in  Figures  6 and  7.  As  expected,  ■ -M,,  and 
agreement  is  quite  good  between  results  in  Fig.  6 and  those  obtained 
for  d » 0 (8).  Gain  and  normalized  far  field  patterns  for  these  cases 
are  given  in  Figures  8 and  9. 

Two  examples  done  by  Neerhoff  and  Mur  [4)  for  a slit  with  w ' lu 
(micron  ),  d ■ .lp,  and  A •*  .4353p  appear  in  Figures  10  and  11  when 
region  b contains  free  space  and  regions  a and  c have  different  per- 
mittivities. The  magnitude  of  currents  M0  is  compared  and  agreement 
is  quite  good.  Figures  12  through  15  show  the  effects  of  having  a 

lossy  medium  in  region  b for  a slit  with  w - 1\  , d - .25A  , when 

a a 

regions  a and  c are  free  space. 

Figure  16  shows  gain  and  normalized  far  field  patterns  for  a 
slit  with  varying  thickness  and  different  values  of  e^.  Our  results 
agree  well  with  the  same  example  computed  in  [4).  This  case  was  also 
experimentally  measured  in  [9|  with  slight  discrepancies  in  the  mag- 
nitudes of  the  sidelobes  and  nulls  when  compared  with  our  results. 
Figures  17  through  19  show  the  magnitudes  and  phases  of  magnetic  cur- 
rents computed  for  the  slits  in  Fig.  16. 


Figures  20  and  21  show  the  effect  on  transmission  coefficient 


for  various  slits  when  the  source  is  placed  at  different  angles  of 

incidence  as  measured  from  the  negative  x-axis.  This  transmission 

coefficient  is  different  from  Eq.  (44)  in  that  P,  is  now  the  time 

inc 

average  power  intercepted  by  the  slit  when  the  source  is  at  normal 
incidence  which,  for  plane  wave  simulation  as  we  have  done,  amounts 
to  multiplying  Eq.  (44)  by  cos<J>  where  4>  is  the  angle  of  incidence. 
This  is  done  to  facilitate  comparisons  with  results  given  in  [12] 
and  [13]. 


POSITION 


Fig.  7.  Magnitude  and  phase  of  Mj  (squares)  and  M 

for  slit  w » .4\  , d ■ .001\  , k,  ■ k “ 
a a b a 

- 3u  ; c)  u - 10y  . N 


2 (triangles) 

k and 

o 

10. 


a)  u 


U : b)  u 


-90  -60  -30  0 30  60  90 

ANGLE 

Fig.  8.  Gain  and  normalized  field  patterns  for  slits  in  Fig.  6. 
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MAGNITUDE 


POSITION 


Fig.  14.  Magnitude  and  phase  of  (squares)  and  M, 

fore  ■ C ■ C » d ■ . 25X  , w • IX.  and 
a c o a a 

r - (1-J10)c  . N - 20. 


( t r langlos) 


TERN 


MAGNITUDE  PHASE 


■ 


Fig.  17.  Magnitude  and  phase  of  (squares)  and  (triangles) 
for  slit  a in  Fig.  16.  N ■ 30. 
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ANGLE  OF  INCIDENCE 


20.  Transmission  coefficient  times  cos<}>  versus  <f>  where  <J>  is 
the  angle  of  incidence  measured  from  the  negative  x axis 
for  w ■ .81  . d ■ .25 


1.0 


0.5 


0.0 


ANGLE  OF  INCIDENCE 


21.  Transmission  coefficient  times  cost))  versus  <}>  where  $ i8 
the  angle  of  incidence  measured  from  the  negative  x axis 
for  k ■ k,  ■ k ; a)  w ■ 1.02X  ,d  • 10  X , b)  w ■ .511  , 
d - l6"AX  7 Squires  and  triangles  represent  values  taken 
from  [13]a  for  d ■ 0. 
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VIII.  DISCUSSION 


A general  purpose  computer  program  has  been  developed  for  calcu- 
lating the  transmission  characteristics  for  a slit  In  a conductor  of 
finite  thickness  separating  two  different  media.  The  slit  may  be  filled 
with  lossy  media  as  specified  by  the  user.  It  is  shown  In  Appendix  A 
that  this  formulation  reduces  to  the  case  of  a slit  in  a conducting 
screen  as  the  conductor  thickness  d approaches  zero.  However,  when 
media  b is  lossless  and  the  slit  dimensions  w or  d are  such  that 


( V + (7) 2 - (7-) 2 
w d Xb 


is  satisfied  for  any  integer  p and  m,  some  terms  of  [Y  ] become  infinite 
and  our  solution  falls.  In  this  case  we  obtain  a solution  by  perturbing 
the  dimensions  w or  d slightly  so  that  the  region  b is  no  longer  a resonant 
cavity. 

When  using  the  computer  program  to  check  with  results  given  for 


slits  when  d ■ 0,  a relatively  small  -slue  ot  d may  be  used  (i.e.. 


d 10  ).  However,  because  the  matrix  composed  of  [Y^Sa] , [ y'5 ] , and 

fi 


,hsc. 


[Y  ] becomes  more  ill-conditioned  as  d approaches  zero,  results  obtained 
for  extremely  small  values  of  d may  be  questionable. 

The  examples  seem  to  indicate  adequate  convergence  when  the  number 


of  expansion  functions  is  increased.  Also,  in  computing  the  infinite 

,.b 


sums  involved  in  determining  [Y  ] , no  more  than  50  terms  were  used  for 
any  of  the  exeaples,  since  increasing  this  number  did  not  appreciably 
affect  the  results. 
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Appendix  A 


NOTES  ON  THE  COMPUTATION  OF  [Y  ] 


The  elements  of  the  eubmatrlces  of  [yj  ere  given  by  equa- 
tlons  (32) , (33) , (40) , and  (41) . All  of  the  equations  have  the 
following  Integral  In  conson: 


I - | | cos  cos  dy  dy' 


(A-l) 


Aal,2  Anl,2 


where  Am^  ^ An^  ^ are  Intervals  over  which  the  expansion  or  testing 


function  is  nonzero  on  or  for  m,  n-l,2,...,N.  Using  the  trigono- 


itrlc  identity 


sin  pnx  - sin  p(n-l)x  ■ 2 sin  cos  |p  (n  - j)xJ 


where  p and  n are  integers,  Eq.  (A-l)  becomes 

. 2 
sin 


2 £TT 


1 “ Uy)2  C°B  [“N1  <m  ' ¥\  C°8  EnL  (n  " 2}] 


(A-2) 


k2N 


where  Ay  ■ -jf-.  Now  write  matrices  [Y11]  and  [Y12]  as 
N 


IY11]  - l a [F^p;] 


(p) 


*•  p 

p-o 


(A- 3) 


and 

where 


[Y12]  - l b[FVP'] 


(p) 


(A-4) 


P"0 


>2  . . . 2 £ir 


Jue.  c (Ay)  cot  k d sin 
b p XE_ 


_2H 


W if/ 


(A-5) 


Juc.  c (Ay)2  cac  k d sin2 
b P 5E 2M 

k . 2 

xp  Q 


_E_ 

w 


(A-6) 
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The  matrix  [F^]  is  given  by 


[F(p)]  - f ? 

P P 


(A- 7) 


where 


f - [cos  cos  3-^-, cos  2 
P 2N  2N 


pn(N  - j) 


(A-8) 


Matrices  [Y12]  and  [Y22]  may  be  written  in  terms  of  the  above  according 


to  Eq.  (16) . 


It  is  evident  that  for  various  values  of  k or  k d when  the 

xp  xp 


medium  in  region  b is  lossless,  the  coefficients  a and  b become  infinite. 

P P 


These  cases  are  summarized  as 


1)  d - o,  kxp  i 0. 


2)  k d mir  for  m ■ 0,  1,  2,... 

xp 


To  examine  these  cases  in  further  detail  we  rewrite  Eqs.  (11a, b)  using 


Eq.  (16): 


[Y11  + Y*18®]^  + [Y12]V2  - I1 


(A-9) 


[Y12]V1  + [Y11  + Yh8C]V2  - 0 


(A-10) 


For  case  1 when  d is  very  small,  a suitable  approximation  to 


cot  x and  esc  x yields 


[Y1-1  j » [Y12  ] 


[Y11]"1  - [0] 


(A-ll) 


iV.  *. 


Next  multiply  Eq.  (A-10)  by  [Y**]  * to  obtain 

Vx  - - V2  (A-ll) 

than  aubtract  (A-10)  from  (A-9)  and  use  Eq.  (A-ll)  to  get 

[Yhaa  + yhacj-^  _ ji  (A-12) 

This  is  the  expected  result  when  the  thickness  of  the  conducting 
screen  Is  zero  [5,  p.  7]. 

An  analytical  expression  for  Eqs.  (A-9)  snd  (A-10)  as 

k d -*■  mw  in  case  2 Is  not  attempted  here.  Instead,  when  the  dimensions 
*P 

v and  d ara  auch  that  case  2 arises  for  some  integers  m and  p,  the  com- 
putations are  dona  for  dimenalons  w + AX^  or  d + AXfe  wherever  appropriate 
where  AX^  is  a small  fraction  of  the  wavelength  in  region  b. 
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Appendix  B 
COMPUTER  PROGRAM 


B-l.  Required  Data 

The  required  date  cards  are  read  in  according  to  the  format: 

100  FORMAT(6E11.4) 

101  FORMATC6I5) 

READ(l.lOO)  NMUA,  NEPSA 
READ( 1,100)  NMUB,  NEPSB 
READ(l.lOO)  NMUC,  NEPSC 
READ(l.lOO)  W,  D 

READ(1,  100)PHIO,  DPHI , FMC 
READ(1,101)N,  NI,  NT,  NEXC 
READ(l.lOl)  I CUR,  IGA 
DO  1 I - 1,  NEXC 
READ(l.lOO)  XSC(I),  YSC(I) 

1 CONTINUE 

The  input  parameters  are  defined  as 

NMUA  - y /u 

a o 

NEPSA  - e U 
a o 

NEPSB  - e.  U 
b o 

NMUB  - y.  /y 
b o 

NMUC  - y /y 
c o 

NEPSC  - t /e 
c o 

W - slit  width  in  meters. 

D - slit  thickness  d,  in  meters 


44 


PHIO  - angle  In  degrees  at  which  first  gain  measurement 
Is  computed. 

DPHI  - Increment  In  degrees  at  which  gain  pattern  Is 
computed. 

FMC  - frequency  of  magnetic  line  source  In  megahertz. 

N • number  of  expansion  functions  on  slit  face  or  ^ 
(total  number  of  unknown  magnetic  currents  - 2N) . 

NI  - number  of  gain  measurements  to  be  computed  for  each 
pattern. 

NT  ■ number  of  terms  used  to  approximate  the  infinite 
summations  in  Eqs.  (A-3)  and  (A-4) . 

NEXC  - number  of  excitations  (i.e.,  number  of  right  hand  sides 
to  Eqs.  (11a,  b)). 

ICUR  ■ Integer  option  variable.  If  ICUR  - 1,  the  magnetic 
currents  will  be  printed  for  each  excitation.  If 
ICUR  i 1,  printout  will  be  bypassed. 

IGA  - integer  option  variable.  If  IGA  - 1,  a gain  pattern  and 
far  field  patterns  will  be  computed  for  each  excitation. 
If  IGA  J 1,  this  computation  will  be  bypassed. 

XSC(J)  ■ x coordinate  of  jth  magnetic  line  source  in  meters. 

YSC(J)  ■ y coordinate  of  Jth  magnetic  line  source  in  meters. 

The  minimum  storage  allocations  in  the  main  program  are  given 
by: 


i 


L 
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COMPLEX  Y (N*(2*N+1) ) , VM(2*N),  YHSA(N) 
YHS  (N),  VM2 (NEXC*2*N) 


DIMENSION  GA(NI) , FP(NI),  PHI(NI), 

PHIR(NI) , XSC(NEXC),  YSC(NEXC) 

For  subroutine  TRANS  the  minimum  allocation  Is 
COMPLEX  YAUX(N) 
for  subroutine  GELS, 

COMPLEX  P.(N*2*NEXC)  , AUX(N*2*NEXC) 
and  for  subroutine  YB, 

COMPLEX  STl(NT),  ST2(NT),  CXP(NT) 

DIMENSION  SINC2 (NT) 
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B-2.  Main  Program  Description 


i 

Li 


The  first  quantities  to  be  computed  are  the  wavenumbers  and 
Impedances  In  each  region  normalized  by  kQ  and  nQ  respectively  where 
kQ  ■ wavenumber  of  free  space  and  nQ  “ impedance  of  free  space.  It 
has  been  noted  that  whenever  the  dimensions  w or  d are  integral  multiples 
of  a half  wavelength  In  region  b (for  the  lossless  case)  computational 
difficulties  are  encountered.  This  condition  is  checked  for  in  statement 
50  where  the  dimensions  are  perturbed  slightly  if  necessary.  Next  the 
dimensions  are  changed  to  electrical  lengths  by  multiplying  by  kQ.  The 
actual  ordering  of  the  magnetic  current  expansion  functions  along  slit 
faces  and  T 2 shown  in  Fig.  B-l.  Here,  (Ay^,  Ay2,...,  Ay^}  are  the 

regions  on  which  + 0 for  n-1,  2 N and  AyN+2*  *'*’  Ay2N^ 

are  the  regions  on  which  M2  ^ 0,  for  m - 1,2,...,N. 

The  first  step  in  solving  Eqs.  (11a, b)  is  to  compute  the  necessary 

matrix  elements  as  given  by  Eqs.  (12).  For  convenience,  each  side  of 

Eqs.  (11a, b)  are  multiplied  by  the  factor  r^k^.  Since  [Y^sa]  and  [Y*180] 

11  12 

are  symmetric  Toeplltz  matrices  and  [Y  ] and  [Y  ] are  symmetric,  only 
the  upper  right  triangular  portion  of 


[Y]  - n k 
o o 


[Y11  + Yhsa]  [Y12] 


[Y12]  [Y11  + YhaC] 


(B-l) 


is  computed.  This  is  stored  in  array  Y in  the  main  program  by  columns. 

The  first  column  of  n k [Y^sa]  is  computed  from  Eq.  (18)  with  statement 

o o 

51  and  the  result  is  stored  in  array  YHSA.  The  first  column  of 


_ I 
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hsc 

r|  k [Y  ] is  next  computed  and  stored  in  array  YHSC  in  statement  52. 
o o 

Statement  53  computes  the  upper  right  triangle  of  the  matrix 


and  stores  the  result  temporarily  in  array  Y.  DO  loop  2 then  adds 
hsd  hsc 

r|  k [Y  ] and  r|  k [Y  1 to  the  proper  elements  of  the  above 

o o o o r r 

temporary  matrix  and  stores  the  final  result  in  array  Y which  is 
the  upper  right  triangular  portion  of  Eq.  (B-l). 

Each  excitation  vector  n corresponding  to  a magnetic 

line  source  at  (XSC(I),  YSC(I))  is  computed  in  statement  54  and 
stored  temporarily  in  array  VM.  This  operation  is  done  NEXC  times 
in  DO  loop  3 where  each  VM  is  stored  consecutively  in  array  VM2. 

The  solution  to  Eqs.  (11a  and  b)  is  then  found  by  statement  55  where 
the  input  array  VM2  now  contains  the  solution  for  each  excitation. 

The  coefficients  of  magnetic  currents  are  printed  out  in  DO  loop  4 
of  ICUR  “ 1.  DO  loop  17  computes  a transmission  coefficient  for  each 
excitation  and  a gain  and  normalized  field  pattern  if  IGA  * 1. 
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f 


n r\  r> 


MAIN  PROGRAM  LISTING 

THIS  PACE  IS  BEST  QUALITY  PRACTICABLE 
KROM  COPY  FURtUSh£D  TO  ilDC 


SJOB  AUCK .T IME=2.RAGES=40 

c CCRRECTEO  TRANSVERSE  ELECTRIC  THICK  SLIT  PROGRAM 

C EASED  ON  REPORT  TIR-77-9 


COMPLEX  VM2<  1820)  «Y<  528)  ,YHSA<  16)  .YHSC<  16)  «VM(i32) 

COMPLEX  CSQRT .CSI N.CCOS. HANK02.HANK1 2.CONJG 
COMPLEX  E T ANB • U • NMU6 . NEPSB , WA  VNB. STK 

DI MENS  ION  XSC (91 ) . YSC(  91 ) • PHI <9  I)  .PHIRI91)  . G A C9 1 ) . FP ( 91 ) 

COMMON  W AVNO. PI . U 

REAL  NMUA.NMUC.NEPSA.NEPSC 

DATA  ETA/376. 730/ . EP S 0/8 . 85E- 1 2/ 

PI  = 2. 141 593 
U=  < 0 . . 1 . ) 

100  FORMATloEll .4) 

101  FORM  AT (615) 

RE  AC ( 1*1 0 0)  NNUA.NERSA 
RE  AC ( 1.100)  NMUB. NEPSB 
REACH. 100)  NMUC.NERSC 
R E A C ( 1.1  00)  W.D 
REACH. 100)  PHIO.CPHI.FMC 
R E AC  < 1.101)  N.NI.NT.NEXC 
REACH. 101)  ICUR.IGA 
DO  1 1=1  • NEXC 
REACH, 100)  XSC(  1 ).YSC(I  ) 

PR  I NT . I . XSC( I ) • Y SC (T  ) 

1 CONTINUE 

PRINT. *NEPSB=* .NEPSB .'NMUB*' .NMUB 
PR  I NT. *NMUA=*  . NNUA  »*NEPSA= • . NEPSA 
PR  I NT • 'NMUC=*  • NMUC . * NEPSC= • • NEP  SC 
PR INT • • W = • . W. • D=  • , D 

PR  I NT  . *Ph  10='  • PHI  0. r0PHI= • .DPHI . • FMC=* .FMC 
PR INT  »*N=,.N» • NI  = • .Nl. *NT=* . NT. *NEXC» • .NEXC 
PRINT, *ICUR=  * , I CUR. ' I GA=  • , I GA 

C COMPUTE  WAV  ENUMBER  OF  FREE  SPACE. 

WAVN0=PI 4FMC/1 50. 

COMPUTE  WAVE  NCSL  OF  REGIONS  A . B . AND  C NORMALIZED 

BY  I*  A V NO  . 

WAVNA=SQRT( NMUA4N  EPS A ) 

WAVNE=CSORT( NMUBANERSB) 

WAVNC=SQRT{ NMUC*NEPSC) 

COMPUTE  INTRINSIC  IMPEDANCES  OF  REGIONS  A.  B.  AND  C 

NORMALIZED  BY  ETA. 

E T A NA  = SQ R T ( NMUA/NEP5  A) 

ETANE  = CSQPT<  NMUB/NEPSB) 

ETANC=SQRT( NMUC/NEPSC) 

PR  INT • WA  VNB • ET  ANB • WA VNO 
M=N*(N+1 )/2 
N2=2*N 

M2*N24<N2H  )/2 
50  CALL  PTB ( W.D.WAVNB) 


SO 


JL 


: 


mis  PA®  IS  SKI  «utttn  PMCHCtfU 

UK*  COP*  FURNISHED  TO  DDC 


C CHANGE  DIMENSIONS  TO  FREE  SPACE  ELECTRICAL  LENGTHS. 

W-kTWAVNO 
D- CA  WAVNO 
DY=  »/N 

C CCMP.  FIRST  RUk  OF  ADMITTANCE  MATRICES  FOR  REGIONS 

C A AND  B...  STCRE  IN  ARRAYS  Y HS  A AND  YHSC  RE  SP . 

El  CALL  Y HS  < W A VN  AADY.N.Y  HSA . WAVNA. ET  ANA) 

52  CALL  YHS  ( WAVNCADY  .N.  YHSC.  MAVNC.  ET  ANC  ) 

53  CALL  Y8< Y .M2. WAVNP.FTANH . W .DY  .NT. C .N > 

DO  2 0 1=  1 .N 

PRINT.YHSAI  I)  .YhSCII  ) 

20  CONTINUE 

C ACD  HALF  SPACE  ADMITTANCE  MATRICES  TO  Y1I  AND  Y22 

C OF  TOTAL  ADMITTANCE  MATRIX  Y... 


K=  1 

DO  2 J=1  ,N 

DO  2 1 = 1 . J 

Y(K)  = Y(K)-»YHSA(J-IM  ) 

Y ( NAJANaN)=Y(  MA  J A NAK  ItYHSCI  J - I ♦ 1 ) 

K = K A l 


2 CONTINUE 

DO  2 1 1=  1 .M2 

PR  I NT  . I . Y I I ) 

21  CONTINUE 

c CCMP.  EXCIIATIQN  MATRIX. 

DO  2 1=1.  NtXC 

X S=  X S C I I ) A WAV  NO 
YS=  V SCI  I ) AWAMNO 


54  CALL  TEEXCI WAVNA.  WAVNAAXS. WAVNA AY  S.N. WAVNA ADY  .VM.N2.ETANA.  STK  . W 
1 AVN  A A W ) 


DO  2 J= l ,N2 

VM2  I J A ( I - I ) AN2  ) = VM( J ) 

3 CCNT  INUE 

C SCLUT  ION  OF  MAGNETIC  CURRENTS. 


55  CALL  GELS(VM2.Y.N2.NEXC.M2> 

IF ( ICUP.NE. 1 I GO  TO  5 

102  FORMAT  I I 5 .4E1 7.7 ) 

103  FORM  AT ( • I * . 15X . • M AGNET IC  CURRENTS  FOR  XS  =*.F8.3,3X. 
A* ANC  YS  =•  . F 8 • 3 ) 

DO  4 J=  1 . NEXC 
WRITEI2.1C3)  XSCI J). YSCI J ) 

DO  4 1=1. N2 

K= I ♦ ( J-l ) AN2 
AI  = A IM AG ( VM2<  K)  ) 

AR=FEAL( VM2IK)  ) 

PHASE=180.AATAN2(AI.AR)/P1 
VA  = CABSI VM2IK)  ) 

WRITEI3.102)  I . VM2IK ) . V A. PHASE 

4 C 0 N T I N Ut 

5 CONTINUE 
DOC  I* 1 ,NI 
PHI(II=PHIOA(  t-l  I A DP  MI 
PHIR  < I lap  I A PH  I ( 1 ) /I  80. 

5 CONTINUE 


SI 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
raOil  OOPY  FURNISHED  lOUDC  


C CONP.  OF  TRANSMISSION  COEFFICIENT. 

10S  FORMAT!* I ' . 18X , ' XS ' • 1 8X. • Y S*  t 1X»* TRANSMISSION  COEFF.*) 

WRITE! 3. l CS) 

DO  7 J=1.NEXC 

DC  € 1=1, N2 

VM!  I)  = VM2!I4! J-l )*N2) 

8 CONTINUE 
XS*XSC! J )*HAVNO 
YS*YSC! J) 4WAVN0 

56  CALL  TRANSIN, N2 . Y H SA .PT , T . VM. STK. W. X S , YS .ETA NA  . WAVNA . TN  ) 

108  FORMAT!*  *.4E20.7) 

WR I TE !3  » 1081  XSC ! J )• YSC! J 1 , T , TN 
MR  I T E ! 2 • 1 CO ) XSC! J) .YSC! J) ,T,TN 
IF!  ICA.NE.l > GO  TO  9 
C --COMP.  CAIN  PATTERN. 

57  CALL  GAIN1PHIR.NI  .N-.N2.VM.DY.W.GA  .WAVNC.ETANC.PT.FP) 

ICC  FORMAT ! F 2 0 • 1 .2E20.71 

107  FORMAT! • 1 *, ‘POWER  GAIN  PATTERN  IN  REGION  C ANO  ANGLE  IN  DEGREES') 
WRITE (3. 1 07) 

DC  10  1=1. N1 

WRITE13.1C6)  PHI! I).GA! I ) ,FP! I) 

10  CONTINUE 

9 CONTINUE 

7 CONTINUE 

STOP 

END 


SAMPLE  INPUT-OUTPUT 


SCATA 

1 -0.10  000  0 OE  03  0.4000000E  00 

NEPS6=  i 0.100000CE  01  , O.OOOOOOCE  00) 

NMUB=  ! 0.1000000E  01  , O.OOOOOOOE  00) 

NWU A=  0. 100000CE  01  NEPSA=  0*  1000000E  01 

NMUC=  O.ICOOOOOE  01  NEPSC=  OrflOOOOOOE  01 

W=  0.8000000E  OC  C=  0.2500000E  00 

PH I 0=  -C.9000000E  02  CPHl*  O.IOOOOOOE  02 

N*  10  N I = 19  NT=  50 

ICUR=  1 ICA=  1 

FMC=  0.3000 0 OOE  03 

N E XC=  1 

! C.100000CE  Cl  • C.OOCOOOOE  00) 

! O.IOOOOOOE  01  . O.OOOOOOOE  00)  0.6283185E  01 
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THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 

PROM  COPY  FURNISHED  TO  DDC • 

PRlNTTUT  or  FIRST  column  or  YmSA  AND  Yf«C... 


0.  1250  8 206  00 

• 

0.16588S5F  00> 

< 0*123082 

02 

00  . 0.1  63.3633  f 00 

1 

O.UT?»7*r  00 

• 

0.645736  62-01  I 

1 0 . 1 1 7?  7 7*r 

oo  . o.6*373eer-oi 

1 

0.4591 6 062-01 

9 

-0. 024086 4F-02 1 

< 0.0551 5062 

-0  1 . - 0.47008842—0? 

1 

0.6  373  6 649-01 

• 

-0.45075742-01  1 

( 0.6373  5842 

-01  . -0.46075742-01 

) 

0.27641 902-01 

• 

-0.2  3124842-01 1 

< 0.27601002 

-01  . -0.631746*2-01 

) 

— 0.62745692-02 

• 

-0.61317362-011 

< -0.62705652 

-02  . -0.61317  36F  — 01 

1 

— 0.325? 630 F -01 

9 

— 0.4  52  530  8F— 01  1 

1 -0.32536303 

-01  . -0.4 67 53 062  — 01 

) 

-0.4T22714F-0I 

• 

-0.22742572-01 1 

1 -0.472271*2 

-01  . — 0.  7 2 74 2 57 F— 01 

) 

-0.4548727F-01 

• 

0.28703642-021 

< -0.48487J72 

-01  . 0.787036*2-02 

1 

-0.39093 I7F-01 

• 

0 .24702462—01  1 

1 -0.3909317c 

-01  . O.747O746F-0I 

) 

PRINTOUT  nr  UPPER  PT.  TRIANGLE  OF 

V.  .. 

t 

0.19605202  00 

t 

0.38441302  001 

58 

< 

0.00000002 

00 

9 -0 

.I53I661F  00) 

2 

0. 11T2TT4F  00 

9 

0 .70407342—0 1 1 

30 

< 

0.00000002 

00 

9 -0 

.1607446?  90) 

3 

0.12508202  00 

• 

0.2746003*  001 

60 

< 

0.00000002 

00 

9 -o 

.153*1*7?  00) 

A 

0.9951 698F -01 

• 

—0.11 207 I 7F  001 

61 

( 

0. OOOOOOOF 

00 

9 -0 

.1040408?  00) 

S 

0.11727742  00 

• 

0*25*44012-01  1 

62 

( 

0. OOOOOOOF 

00 

9 -o 

.3081 1 °6F -01 ) 

3 

0.1 2506202  00 

9 

0.27665302  001 

63 

C 

0. 00000002 

00 

9 0 

.402  43  50? - 01  ) 

T 

0.6373584F-01 

• 

-0.10560032  00 1 

64 

1 

0. OOOOOOOF 

00 

9 0 

.11677802  00) 

9 

0.95516482-01 

9 

-0.10012002  001 

65 

( 

0.00000002 

00 

9 0 

•1557757F  00) 

9 

0.11 727 742  00 

9 

0.60006802-01 1 

65 

< 

0.12506202 

00 

9 0 

.38*4130*  00) 

10 

0.12506202  00 

• 

0.34078202  001 

67 

1 

0.00000002 

00 

9 -0 

•11 48069F  001 

11 

0.27641  OOF-01 

9 

-0.20741515  001 

66 

( 

0. OOOOOOOE 

00 

9 -o 

.1766464*  001 

12 

0.63735442-01 

9 

-0.15124642  00) 

60 

( 

0.0000000* 

00 

9 -0 

. 1334733F  00) 

13 

0.95516982-0 1 

9 

-0.37407302-01 1 

70 

c 

0. OOOOOOOF 

00 

9 -0 

. 13  7338  7F  00) 

1* 

0.1 1727742  00 

9 

0.14008652  OOl 

71 

< 

0.00000002 

00 

9 -0 

.1106703?  00) 

IS 

0.12508202  00 

9 

0. *1646175  00) 

72 

( 

0. OOOOOOOF 

00 

9 -0 

.601 8623F-0 1 1 

13 

-0 .6  2795655  — 02 

9 

-0.15144502  OOl 

73 

< 

0.00000002 

00 

9 -0 

• 234  657  7* - 01  1 

IT 

0.2  769 1 90F— 0 1 

9 

-0 .1 36784 72  001 

7* 

« 

0*00000002 

00 

9 0 

.367738 1F-01 ) 

19 

0.6373584E-01 

9 

-0.72066372-011 

75 

( 

0.  OOOOOOOE 

00 

9 0 

. 877301 7F -01 ) 

19 

0.98516482-0 l 

9 

0.20181502-011 

75 

( 

0 .OOOOOOOE 

00 

9 0 

.1167780?  001 

>0 

0.11727742  00 

9 

0.18745632  00) 

77 

1 

0.  1 177774E 

00 

9 0 

.7  0*  47 3 OF -01) 

21 

0.12506202  00 

9 

0.416  061*42  00> 

78 

< 

0. 1250620F 

00 

9 0 

.77*80432  00) 

22 

— 0.32536302— 01 

9 

-0.75*50012-01 1 

70 

< 

0. OOOOOOOF 

00 

9 -0 

.15316812  00) 

23 

-0.62793652-02 

9 

-0.82764332-01  1 

80 

( 

0. OOOOOOOE 

00 

9 -0 

. 1 33473JE  00) 

24 

0.2'r69l  902-01 

9 

-0.60603632-011 

61 

c 

0.00000002 

00 

9 -0 

.11061712  00) 

29 

0.63735842-01 

9 

-0  .33607*12-01 1 

62 

1 

0. OOOOOOOF 

00 

9 -0 

.665P09OF-0 1 1 

26 

0.95516982-01 

9 

0.201618 32 -Oil 

83 

< 

0*00000002 

00 

9 -0 
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PRINTOUT  3 F EXCITATION  MATRIX... 
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-0. 10062 99E 

0 1 
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-0.1 0053  04p 

0 1 

0 . 1 76  07  2 2F-0  2 ) 

3 ( 

-0.1 0O5307F 

0 1 

0.61 3592 1 E- 03 ) 

4 ( 

-0.1005307E 

0 1 

-0 .1653B22E-03) 

5 < 

-0. 1 0053 07F 

0 1 

-0.61 3592 1 E-03) 

r>  ( 

-0.1 00530TF 

0 1 

-0.61 35921F-031 

7 < 

-0 • 1 005307F 

0 1 

-0 . 1 653822F-031 

3 ( 

-0. 1 005307F 

01 

0.61 35921 E-03) 

9 < 

-0.1 005304F 

01 

0.1 76  07j  2E - 02 1 
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0. OOOOOOOE 
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70  ( 

0. OOOOOOOF 
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0. OOOOOOOF  00) 

MAGNETIC 

CURRENTS  FOR 

1 XS 

“-10O.000 

4Nn 

YS  = 0.400 
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01 

-0.1 21 4056F 
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00 
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00 

- 0. 1 75  853OF 
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-0.9676579F 
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00 
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1 4 
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00 
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1 5 
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02 

-0  .91  B376‘*F 

00 
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00 

-0.904561  pc 

07 

1 6 

-0. 731 2000F- 

02 

— 0.9  1 8367  OF 

00 

0.91 8396 1 F 

00 

-0.904561  OF 

0? 

17 

-0*  21 56668F- 

01 

-0.9213762E 

00 

0. 9216785“ 

00 

-0.  91  3407TF 

07 

i a 

— 0 • 59  3 1 2 4 0 F - 

0 1 

-0  .9  31  1371  F 

00 

0 • 9 3 30  7 43F 

00 

— 0.  9364  4 71  f 

0? 
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-0.9  136 149F 
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xp  YS  TRANSMISSION  OFFF. 

-0.  lOOOCOCE  03  0.4  0000  OOp  00  0.9396070F  00 
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POWER  GAIN  PATTERN  IN 

REGION  C AND  ANGLE 

IN  DEGREES 

-90  .0 

0. 1 042085F 

00 

0. 51 0009^F-01 

-80.0 

0.  1200446E 

00 

0.54  73  917  E— 0 1 

-70  .0 

0.  176271 7F 

00 

0*  66  33 1 2 1 F- 01 

-6  0.0 

0.  29  8677  1 E 

00 

0.86  3431  IE  — 01 

-50.0 

0. 5253934E 

00 

0.  1145167F  00 

-40.0 

0. 8892853E 

00 

0.14  89866E  00 

-30.0 

0. 1383474P 

01 

0.18 58285 E 00 

-20.0 

0. 1925999E 

01 

0.2192577E  00 

-10.0 

0.2362068F 

01 

0.2428136F  00 

0.0 

0. 400632  5F 

02 

0. 1 OOOOOOE  01 

10.0 

0.2362071E 

01 

0.2428138“  00 

20.0 

0. 1925999E 

01 

0.21 92577F  00 

30.0 

0.  1383476E 

01 

0.185P286E  00 

40.0 

0. 889287  OE 

00 

0.148P868F  00 

50.0 

0. 525393  8E 

00 

0.1145158E  00 

60.0 

0.  2986779E 

00 

0. 8634323F— 01 

70.0 

0. 17  6272  0^ 

00 

0. 6633127p-01 

80.0 

0.  12  0044  8E 

00 

0.5473922  E-0 1 

90.0 

0. 1042085E 

00 

0.  51  00096-F-01 

COPE  USAGE 
COMPILE  timE= 


OBJECT  CODE=  49704  BYTES*  ARRAY  A9FA  = 45516 

3*56  SEC  «EXECUTI  Qfsj  tjmE:=  4*88  SFC*  WATFIV 


56 


F 


; 


! 

! 

' 

I 


• _ • ■ * 1 a. 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
raOli  COPY  FURNISHED  TO  DDC  — 


B-3.  Descriptions  and  Listings  of  Subroutines 


SUBROUTINE  PTfcM  w.D.WAVNB ) 

CCMMCN  WA VNO.PI .U 

COMPLEX  U • WAVNB 

WLE=2 .*P I/t REAL! W AVNE)** AVNO ) 

DD=.001*WLB 
1 = 0 

IF(O.EG.O.»  GO  TO  IT) 

IF(AIMAG( WAVNB ) • NE  .0  . ) RETURN 
00  1 1=1.20 

IF!ABS!MLE*I/2.-0).LE.l.E-04)  GO  TO  10 

1 CONTINUE 
GO  TC  11 

10  0=  C *C 0 

100  FORMAT!* -*,*D=  * . 1 5.3X  . • T I ME  S HALF  M A VELENiGT  H IN  REGION  B.  CHANGED 

* TC  * . E 14 . 7 ) 

NR  I T E ( 3 . 1 00  ) 1.0 

11  IF!AIMAGi  feAVNB  ) .NE  .O  • ) RETURN 

00  2 1 = 1 .20 

IF ( AES(WLE* 1/2 . ) .LE*  l.E-04  ) GO  TO  15 

2 CONTINUE 
GO  TO  16 

15  w=M-»DD 

150  FORMAT!*-*. *W=  * . I5.3X. *T I MES  HALF  WAVELENGTH  IN  REGION  B.  CHANGEO 

* TO*  .E14.7) 

WR  I T E ! 3 . 1 50  ) I.W 

16  RETURN 
END 


This  subroutine  checks  to  see  whether  w or  d ere  Integral 
multiples  (up  to  20)  of  a half  wavelength  In  region  b.  If  either  d 
or  w satisfy  this  condition  they  are  changed  by 

d - d + .001  Ab 
w - w + .001  A. 

D 

If  region  b Is  lossy,  the  parameters  d and  w are  unchanged  unless 
d ■ 0 In  which  case  only  d Is  changed. 
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C LISTING  CF  SUBROUTINE  YB 

C 

SUBROUTINE  YB < Y . M2 .W AV NB . ETANB • M. CY  . NT . D , N ) 

COMMON  In  A VNO  • PI  • U 

COMPLEX  Y ( M2 ) • YB 1 1 .YBl 2. YSUM 1 I , YSUM1 2 .KBD. U, CONST. WAVNB. ETANB 
COMFLEX  ST I ( 150) • ST2< 150) .CXP( I 50  ) , C SORT ,C S I N . CCOS 
DIMENSION  S I NC2( 150) 

N 2=  N*  2 
M=N* ( N4l )/2 

CO  N£T=-<U*OY*DY/W/ ETANB 
KBC  = VAVNB4D 

IF (ABSIAI MAGI  KBO)  ).GE.50.)  GO  TO  I 
YBl  2 = CON  ST/CSINIKBCI 
YB  1 1 = Y B 1 2 4CC0  S ( KB  0 ) 

GO  TC  2 

1 MF=-1 

IF< A IMAG< KBD) . GT. 0. ) MF= 1 
YB 1 1=MF*CCNST*U 
YBl  2=0. 

2 CONTINUE 

DO  2 1=1.  NT 

CXP  ( I )=C  S ORT ( 1 I*PI/W/WAVNB)**2  ) 

SINC2(I)=(SIN(  I*  F 1/N2 ) /(  I *P I /N 2 ) ) 4* 2 
K B C=  C XP ( I )*MAVNB*D 

IF ( AES < AIMAGI KBD ) ) .GE.50. ) GO  TO  4 

ST  111  )=CCCS(KBD)*SINC2(I ) /CS INC  KBC )/CXPC  I) 

GO  TO  e 
4 MF= - 1 

IF ( A IMAGI KBD) .GT.O.)  MF= 1 
ST  1 ( I ) =S I NC2 ( I )/CXPC  I)*MF*U 

6 CONTINUE 

3 CONTINUE 

DO  7 1=1, NT 

KBO=*AVNB*CXP< I )*D 

IF<AES<AIMAG< KBO) ).GE.50.)  GO  TO  € 

ST 2 < I ) =ST  1 < I)/CC0S<KBD) 

GO  TO  9 

8 ST  2 < I ) = 0 * 

9 CONTINUE 

7 CONTINUE 
K = 1 

DO  10  1= 1 *N 
DO  10  J=  1 .1 
YSUM1  1=0  • 

Y SUM 12=0 . 

DO  11  I P * 1 , NT 

CS=COS< IP*PI*< I-.5)/N)*COS< IP*PI*< J- .5)/N) 

YSUM1  1=Y S UM 1 1 + ST  1 <IP)*CS 
YSUM 12=Y  SUM  1 24 ST 2 < IP )*CS 
11  CONTINUE 

Y < K ) = YB 1 14YSUM1 l*CONST*2. 

Y<N4< 1-1 ) 4N4 K ) =YB 1 2+YSUM12*CONST*  2 . 

K=K4  1 

10  CONTINUE 


/ 
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K=N_!  from  COPY  FURNISHED  TO  DDC 

N 1 = K 

00  14  1=  1 ,N1 

00  1£  J= 1 »k 

Y(M*(I-l  ) 4N  + J ♦ I * { 1*1  )/2)  = Y(M+(j  + I— 1)*N+(J+I— 1)*(J+I)  / 24  I ) 

15  COM  I NOE 
K=K  - 1 

14  CONTINUE 

K=  1 

OC  1C  1= 1 ,N 
DO  lfc  J=  1 , I 
Y ( mi4N*K  ) = Y(K) 

K-K  •*  1 

16  CONTINUE 
RETUPN 
ENC 


This  subroutine  computes  the  upper  right  triangle  of 


n k 

o o 


and  returns  the  result  in  array  Y.  The  major  task  is  to  compute 


[YU] 

[Y12] 

[Y12] 

[ Y 1 1 ] 

approximations  to  the  infinite  sums  in  Eqs.  (32)  and  (33).  Re- 
writing Y^  as  it  appears  in  Eq.  t32),  we  have 
mn 


. V11 
n k Y 
o o mn 


jk  n (Ay)  cot  k.d 
o o b 


j2ko no(Ay)2-  C°t\J1  " (kbw}  V sin2(g) 


o o_ 

WTL 


p=l 


1 - (t^1-)^ 
k.  w 
b 


(£21)2 


COS 


R? (m  - ¥\ 


cos[f  (n  - |)] 


(B-2) 


where  Eq.  (A-2)  has  been  used.  A similar  expression  is  obtained  for 

12  2 
H k Y . The  complex  constant  CONST  is  set  equal  to  -i(Ay)  k n /wn,  . 
o o mn  J 7 o o b 

Since  k^  is  a complex  number,  we  write  k^  = k^  - jk^  and  if  k^d  50, 

then  the  limit  of  cot  k^d  as  k^d  ■*  °°  is  used  in  computing  the  p = 0 term, 

YBll,  for  n k Y^.  Next,  DO  loop  3 computes  NT  terms  of  the  infinite 
o o mn 

summation  for  Y^  where 
mn 
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used.  PO  loop  7 computes  the  terms  of  summation  for  n k Y nnd  stores 

o o mn 


them  In  array  ST2  as  ST2(p)  ■ ST1 (p) /cos(CXP (p)k^d  ) . DO  loop  10  computes 
NT  terms  of  Eq . (B-2)  for  m,n  - 1,2,...,N  and  stores  the  result  In  the 
first  N(N+l)/2  locations  of  Y.  Note  the  result  Is  an  upper  right  tri- 
angular matrix  stored  columnwise.  The  upper  right  triangle  of  the 
12 

matrix  n k Y is  also  computed  here  and  stored  in  the  appropriate  loca- 
o o mn 

12 

t ions  in  Y.  Next,  since  ti  k^(Y  ] is  symmetric,  DO  loops  14  nnd  15  fill 

12 

in  the  rest  of  that  portion  of  Y which  n^k^lY  1 occupies.  Finally, 
since  [YU  ] - [Y“],  DO  loop  16  fills  in  the  rest  of  Y which  is  then 
returned  to  the  main  program. 


r 

r 


Lir-TING  OF  SUBROUTINE  YHS 


THIS  PASS  IS  REST  QUALITY  PRACTICABLE 
raOM  COPY  FUWUSHHD  TOLDG  - 


SUPROUTINF  YH S< OY . N , YHS A , W AV , ET AN ) 

CPMPl FX  SUM.HANK02 .YHSA(N) .HANK12.U 
COMMON  WAVNO. PI ,U 
OIMFNSION  t<*)  , A ( 4 ) 

DATA  T/- .4305682, . 4305682 I 6ORQ05 •• 1 ftPCpOB/ 

OATA  A/2*. 3478548, 2*«6S21 45?/ 

DATA  NL/3/ 

YHSA  < 1 ) = 0. 

OX  = DY/NL 
nn  i 1 = 1 .? 
no  2 J=1 ,NL 

X1'APS(  ( I -l  )*DY-Dx*  < J-.5)  ) 

X?=X1 *X1 

H0=X1*(  1 ,-X?/0.*( 1 . — X?  ^ 25 « * ( 1 .-X2/4P.  > ) ) 

HleX?/3.*( 1 • - X2/15.*< I ,-X2/35.*< 1 .-X2/ft3. ) ) ) 

Y HS  A ( 1 ) = YHS  A(  1 ) ♦OX  * ( X 1 * ( ( 1 .-HI  )*HANKO?<  XI  ) *HO*HANK  1 2(  X 1 ) ) ) 
? CONTINUE 

1 CONTINUE 

YHSA<  1 )=YHSA< l)/(2 .*FT AN* WAV ) 

C 2 = DY /? . 

DO  3 K=2,N 
C4=(K-.5)*0y 
SOM=Oi 
no  a J= l .4 
00  5 1=1,4 

SUM=SUM*A<  J )*  A<  1 ) 4 HANK  0 2 < ABS  OY  * < T < I ) -T  ( J ) > *C  2-C  A > ) 

5 CONTINUE 

4 continue 

YHSA ( K 1 =SUM*C2**2/ (2.* FT AN* WAV) 

3 CONTINUE 

PFTURN 

f Nr 

This  subroutine  computes  Eq.  (18)  or  (20)  multiplied  by  the 

factor  (n  k ).  The  parameters  are: 
o o 

DY  - Ay  times  k or  k where  Ay  is  defined  in  Fig.  B-l . 

'm  a c m 

listt  hsc 

N - number  of  elements  in  one  column  of  (Y  ] or  (Y  ). 

YHSA  - array  containing  first  column  of  [y*'8a]  or  [Y*1  ]. 

WAV  - k /k  or  k /k 
a o co 

etan  - n /n  or  n /n 

a o co 


ns 

The  exact  elements  stored  in  array  YHSA  are  n k Y ' . 

O O HI  I -I 


i m- l|2|ii«|Ni 


When  m + lt  a four-point  Gaussian  quadrature  formula  Is  used  to  evaluate 

h sa 

the  integral.  Thus,  for  example,  n k Y becomes 

O O HI  | l 


(B-6) 





-•  " 


k Ay,  k Ay  44 


, „ hsa  1 ,V*7lwW  r v . . „(2),,.  . u 

'oV..l  ---  k,  n.  (— )<—)  j,  1 AiAJH. 


k 'xn 
o o 


where  the  A^,  are  weight  coefficients  and  nodes  respectively  and  are 
stored  in  arrays  A and  T.  When  m-1 , the  integrand  is  singular  so  the 
approximation 


. vhsa 
n k Y,_ 
o o 11 


1 NL  f 
I H 

(Vw,J 


2<r><7T>~V*i 

o o 


k Ay  (i  - . 5)  k Ay, 


is  made.  The  integral  inside  the  summation  is  then  computed  with  the 
aid  of  the  identity  [10,  11.1.7] 


H(2)(t)dt  - xH(2)(x)  + ™ (H  (x)  H(2)(x) 
o o i o i 


- H (x)  H^2)(x)} 
1 o 


(B-8) 


where  Hq(x)  and  H^(x)  without  superscripts  are  Struve  functions  of 


orders  zero  and  one  defined  as 


and 


Ho(x) 


H1(x) 


7 


1 .3* 


2 2 2 
1.3.5 


12.3 


2 2 
1.3.5 


2 2 2 

1.3  .5  .7 


- ..•] 


(B-9) 


(B-10) 


C LIFTING  OF  SUBROUTINE  tFFxc 

C 


S'JFPOUT  INF  TFFXC<  WA V.XS ,YS,N, DY , VM, N? .FT  AN. STK ,W  > 
CPPMON  WAVNO.PI.U 

r OVPLFX  VM( N?  ) , STK , SUM, U.HANK02.HANK1 ? 

D1MFNSICN  T(4),A(«) 

DATA  T/- .4305682. .430566?, 1 600906, . 1 609006/ 

DATA  A/2*. 3478548, 2*. 6521452/ 

FS0=S0RT(XS*XS'M  YS-W/2,  )**?> 

S T K = 4 . * U/H A NK 1 2 (RSO)/WAV 
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THIS  PACJE  IS  BEST  QUALITY  PRACTICABLE 
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DY2=Dy/2 • 
no  ? j*i.n 

SUM*  0 • 


SBIS  F&QK  is  BBST  QUALITY  FRACTICABL* 

JHOM  COPT  FURWSHHD  TO  DDO  - 


DO  1 1 = 1 .4 

SUM=SUM*A<  1 )*HANK02  < SORT ( XS* XS+ < ( < J- .*)  + T(  I)  )*pY-YS)  *4  2)  ) 

1 continue 

V M ( J )=-( SUM*STX*DY2 )/( ETAN* 2 • ) 

? CONTINUE 


Nl*  N + 1 

DO  3 I = Nl,N? 
VM( 1 ) =0. 

CONT I NuF 

RFTURN 

FNP 


This  subroutine  computes  the  elements  of  Eq.  (43)  multiplied  by 


kQno.  The  input-output  parameters  are: 


WAV  - k /k 
a o 


XS  - k x 
a s 


YS  - k y 
a-'s 


N • number  of  non-zero  elements  of  excitation  vector. 


DY  - k Ay 
a m 


VM  - output  array  containing  excitation  elements  (right 


hand  side  of  Eqs.  (11)  multiplied  by  n k ) 

o o 


N2  - 2*N 


ETAN  - n /n 
a o 


STK  ■ Strength  of  line  source  which  produces  an  incident  electric 
field  equal  to  unity  at  the  origin. 


W ■ k w 
a 


To  produce  a unit  incident  electric  field  at  the  origin,  the  strength 


of  the  magnetic  line  source  Ku  must  be 

z 


4.1 


STK  - k K - 
o 


/ 2 2 
J x + y 

Vs  ' s 


h{2)  ( /x2  + y2)  x k /k 
1 Vs  1 s sac 


( B— 11) 
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' I 


A four-point  Gaussian  quadrature  rule  is  used  to  approximate  the  integral 
in  Eq.  (43)  thus  the  elements  computed  are 


n k I1 
o o m 


k Ay  n 

- - (-V5)  (j2-)  (k  K) 

i zn  o 
a 


4 

I 

i-i 


ah(2> 

1 o 


/ 


‘V.>2  * 


(V. 


- v2> 


(B-12) 


This  is  done  in  DO  Loop  2. 


c listing  OF  SUBPOUT  inf  TPANS 


SUPPOUTI NE  TP  ANSI N, N2 ,YHSA ,PT ,T, VM, STK ,W,XK ,YK , ET ANA , W AV N A , J N ) 
COMPLFX  YHS A( N> ,YAUX( 30  >»  VM( N2 ) , SUM.S.CCNJG.STK ,U 
COMMON  WAVNO.PI.U 
DO  1 1 = 1 .N 

Y AUX ( I ) =CON JG (YHS  A ( I)  ) 

1 CONTINUE 

Nl=N-l 
SUM=0. 

DO  2 1=1, Ni 
S = 0. 

DO  3 J=i  fI 

Sc  S+Y  AUX  ( I- J*1  ) *CON  JG  ( VM(  N*  J)  ) 

? continue 

NI =N— I 
DO  4 J=1 , NI 

S=  S+Y  AUX  ( J4  1 ) *CON JG  ( V M ( N+  I ♦ J ) ) 

4 CONTINUE 
SUM=SUM4-VM(  N4  I ) *S 

? CONTINUE 

FcO. 

DO  S 1=1 , N 

S=  S + YAUX ( N— I ♦ 1 ) *CON  JG  ( V M( N*  I)  ) 

5 CONTINUF 

S'JM=  SUM+  VM(  N2  )*S 
pTKMe ( C ABS< ST  K ) ) *♦ 2 

A?aW*P 

P2=XK*XK+( YK-W)**? 

C2  = XK*XK  4-YK4YK 

THFT A = A RCOS ( ( B2  + C2  — A2 )/ (2  «*SQRT (B2*C2 ) ) ) 

ThFtAN=2.*ATAN( W/2, /SQRt< XK* XK+YK*YK ) ) 

PT=RE AL (SUM ) 

T=B.*PI*ETANA*PT/( THET A*WAVNA*STKM) 
tN  = 8. I *FTANA4PT/(  Th?tan*WAVNA*STKM) 

FFTUPN 

ENO 
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This  subroutine  computes  the  slit  transmission  coefficient,  T 


according  to  Eq.  (53).  The  input-output  parameters  are 

hisc 

N - number  of  elements  in  first  column  of  [Y  ] 

N2  - 2*N 

tisc 

YHSA  - array  containing  first  column  of  [Y  ]. 

PT  = k n Re  (0  [Yhsc]*  V*} 

O O 2.  I 

T = transmission  coefficient,  Eq.  (53). 

VM  - array  containing  the  solution  to  Eqs.  (11). 

STK  - k K 
o 

W • k w 
o 

XX  “ k x 
o s 

YK  - k y 
o s 

ETANA  - n /n 
a o 

WAVNA  - k /k 
a o 

The  auxilliary  array  YAUX  contains  the  first  column  of  [Y^SC]  . It 

has  a minimum  dimension  of  N.  The  computation  of  Eq.  (50)  is  carried 

out  in  DO  loops  2 through  5 and  the  power  transmitted  k n P , is 

r o o trans 

assigned  to  variable  PT. 


c 

c 


LISTING  OF  SUBPOUTlNF  GAIN 


SUBPOUT INF  GA I N( PH I P » NI  , N, N? • X , DY . W « G A . WA VNC . T TANC « PT * FP ) 
COMPLFx  X(N? ) .U .S . ARG. CFXP 
DIMFNSICN  G A(  91  ) »FP(P1  ) «PHIR(91  ) 

COMMON  WAVNO.PI.U 
0T2*0Y/2. 

DO  2 J>1,NI 
SN*SIN(PHIP(J ) ) 

I F ( SN.FO.O. ) GO  TO  5 
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ft=sinc  Dy;>*WAVNC*SN)/WA  VNC/SN 
GO  TO  6 
FT=1 ./WAvNC 
S=0. 

DO  1 1=1 ,N 

YM=WAVNC*( { I- .5>*DY— W/2. ) 
ARG=U*YM*SN 

S=S+X< I+N)*CEXP(ARG)*FT*2. 

CONT INUE 

F=CABS(S) 

FP( J)=F 

GA  ( J )=F*F*WAVNC/<2  .*ETANC*PT) 
CONTI  NUE 
FM=FP( i ) 

00  3 I=?,Nl 

IF  ( FM.LT.frP(  1 ) ) FM=FP(  I ) 

CONTINUE 

00  A 1=1 ,NI 

FP<  I } =FP(  I ) /FM 

C 0NT I NUF 

FFTURN 


This  subroutine  computes  the  power  gain  pattern  according  to 
Eq.  (64).  The  input-output  parameters  are 

PHIR  = input  array  containing  values  of  <$>  in  radians  at 
which  each  gain  computation  is  made. 

NI  =*  number  of  elements  in  PHIR 
N ■ number  of  elements  in  vector  in  Eq.  (64) 

N2  - 2*N 

X - array  containing  solution  to  Eqs.  (11a, b)  for  a 

particular  excitation. 

DY  = k Ay 
o 

W - k w 
o 


GA  - output  array  containing  G(iJ>) 


WAVNC  - k /k 
c o 

etanc  - n /n 

c o 

PT  - k n Real  {V_  [YhsC  ] * V.*} 
oo  l i 


FP  - output  array  containing  normalized  H-field  pattern. 
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Hi 


*5  ' *«.  • • > 


The  normalized  field  pattern  is 


found  using  Eq.  (63)  and 


FPU) 


H (r  , 40 

m m 


W^ma* 


(B-13) 


C LISTING  of  subroutine  gels 


SUBROUTINE  GEL S C P t * • M, N , M2 ) 

COMPLEX  A(M2) 

COMPLEX  R(  1 820)  .AUX(  2000) 

COMPLEX  P I V I » T0 

00  FORMAT! • 1 • . (WARNING FRROR  CODE  IFR  = *,15) 

EPS=1  .E-07 
IF  <M)24 ,24, 1 

1 IER=0 
P I V=0  • 

L = 0 

DO  3 K=l • M 
L = L+K 

TPA=CABS(A(L) ) 

IF(tFA-PIv)  3 , 3 t 2 
? P I VsTpA 

I=L 
J = K 

3 CONTINUE 
TOL*FPS*PI  V 
LSfsO 
NM=N»-M 
LFND=M-1 
DO  IB  K=1 ,M 
IF  f»I v)24,?4,4 

4 1F(IEP>7.5.7 

5 IF  (PlV-TOL)6,  6,7 

6 I F R = K-1 

7 L T = J — K 
LST=LST+K 
PI  Vl  = l./A(  l ) 

DO  fl  L=K.NM,M 
LL=L+LT 
TP=PI  VI  *R(  LL) 

R ( LL )=R(L ) 

p PCL)=TB 

IF<K-M)9,I9,1  9 

9 LR=LST4<LT* ) )/2 
LL=LR 
L = LST 
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DO  14  II=K,LENO 

L=L+II 

LI  =LL  + 1 

I F(L-LR) l?,  10 ,1  I 
10  A(LL)=A(LST) 

TB  = A(L  ) 

GO  TO  13 
n LL=L+LT 

12  TB=A(LL) 

A(  LL) =A( L) 

1 3 AUX( 1 1 >=TB 

14  A(L)=PIVI*TB 
A<  LST)=LT 

P I V=  0 . 

LLST=LST 

LT=0 

DO  1 B I I =K, LEND 
Pi Vl=-AUX( 1 1 > 

LL=LL*t 

LT=LT+i 

DO  15  LLD=I 1 .LEND 
LL=LL+LLD 

L=LL+LT 

1?  A(L)=A(L)+PWI*A<LL) 

LLSTslLST+I  I 
lr=llst+lt 
TBA=CABS( A( LR ) ) 
IF(TBA-PIV)  17,17,16 
6 PI V=TBA 
I =LP 
J=II+1 

DO  IP  LR=K,NM,M 
LL=LR+LT 

IP  P( LL )=R<LL ) *° IVI*R (LR ) 
ic  I F { LFND  >24,23, 20 

20  1 I =M 

DO  22  1=2, M 

lst=lst-h 
11=11-1 
L=  A ( L ST ) +, 5 
DO  22  J=I1,NM,M 
T B=R ( J ) 

LL  = J 
K=LST 

DO  21  L T = I I .LEND 
LL=LL-M 
K = K+l t 

21  TB=TB- A (K)*P( LL ) 

K=  J+L 

P<  J)=R<<  ) 

22  F ( K) =TB 

23  PETUPN 

24  I ER=  — 1 

WR I TE  <3,100)  IER 
PETUPN 
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This  subroutine  solves  the  set  of  equations 


[A]x  - [R] 

where  [A]  is  a symmetric  matrix.  Input-output  parameters  are 

R = excitation  array  of  dimension  M*N. 

A » symmetric  matrix  of  coefficients,  the  upper  right 
triangular  portions  of  which  is  stored  by  columns. 

M - number  of  equations  in  the  system  to  be  solved. 

N - number  of  right  hand  sides  or  columns  of  R. 

M2  - number  of  elements  in  upper  right  triangle  of  [A]. 

The  algorithm  utilizes  Gaussian  elimination  with  pivoting 
on  the  main  diagonal.  For  further  details,  see  [11]. 


LISTING  CF  FUNCTION  SUBPROGRAM  HA  NIC  1 ? 


COMPLFX  FUNCTirN  H A NK 1 2 ( X ) 

COMPLFX  U/(0..1 

100  fpRMATC  *, ‘WARNING APGUMENT  OF  • . F 1 5 . 4 . 3X  . • HAS  FERN  FNCnuNi-FREP 

1 IN  CDMOUTING  HA  NK 1 2 • / ) 

IFIX.tF.0.)  WPITFI  ? , 1 CO  ) X 
3 SJ 1 = 0. 0 

IF (X .FO.O. ) GO  TC  ? 

7=  ABS(  X) 

IF(7.GT.3.0)  GO  TO  1 
Y=7*7/P.O 

BSJl=x*(0.5*Y*  < -0. 56240985* Y*  <0.21003573* Y*(  -0  • 03054?  8^  *Y  * ( 0.30443 
131 0*Y*(  — 0.00031  761  *Y* 0.00 001109))))  )) 

GO  TO  ? 

1 W=3./7 

FI = . 79  78845  6*  W* I . 00000155  *W *(. 016 50  66 7+w *(.0  9 0171 05*W*< -. 00 24 
195 1 1*W*<  .0011  3653- W*.  00 02 0 0 33  > ) ) ) ) 
ri =.7853981 6*W»  <- . 1 249961 2 *W  * < - . 0 00  056  5 0 + W*  < • 0 067 78 70* W*<  - . OO 
1 0 7434 8* W* < -.0 00795 24+ W*  .00027166))))) 

PS Jl =F 1 *S IN( 2-Pl  ) /SQRT<  Z) 

IF!X.LT.O.O)RSjl =-BSJl 
2 CONTI  NUF 

FSY 1 =_ 1 ,0F75 

IF  ( X.F  O. 0 • ) GO  TO  X 

IFIZ.GT.3.)  GO  tO  4 

THIS  PA 01  IS  BEST  QUALITY  PRACTICABLE 
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RSY1= ( -0 • 6?  661 97  7 4Y* ( 0*2? 120  91 4y*  ( 2. 1 6 827  094Y*  ( - 1 . 3 1 6487  7+  Y*<  0.  ’12 
13951 4 Y*< -0. 04  00 975 4 Y*  0 . 00  27873 )))) ) ) / Z* o . 63 66 1 977 * AL 0G( 0 «5*Z ) *RSJ 1 

GC  TO  3 

4 PSY 1=-F1 *COS( Z-°l ) /SQRT( Z) 

3 CONTINUE 

HANKl 2=RSJl  -U*pSYl 
PFTU°N 
END 

c listing  of  function  suRprogrAm  Hanko? 

c 

COMPLFX  FUNCTION  HANK02(X) 

COMPLEX  U/<  0. .1  . )/ 

13  roBMAT(lH  • ' M AR N I NG  - APGUMF9T  OF • t F 1 * . 4 « 3X . • HA S BFFN  ENCOUNTFPr- 
10  IN  CALCULATING  H ANK02 ( X ) * / ) 

IF<X.LE.0.0)WRITE<3,10)X 

PSJ0=1.0 

IF<X.EO.O.  ) GO  T0  2 
Z=ABS<X) 

IF(Z .GT.3.0130  TO  1 
Y=2*Z/9.0 

PS J0  = 1 ,0+Y* (-2. 24999974Y* ( 1 . 2 5 56 208 +Y * ( - 0 . 3 1 63 P66 +Y *( 0 • 04  4 44  704 Y*  ( 
1-0* 00394 44  4Y* 0.00021) ) ) ) ) 

GO  TO  2 

1 **  = 3.  /Z 

FO  = .7978 84 56*  W*  ( - » 00000  0774**  ( - . 0 055?  74 4 W*  ( - . 0 0 00 951  ? 4 W*  ( . 00 1 
13  72374W* (-. 00 072B054W* .0001  4476 ) ) ) ) ) 

PC= . 7 853981 6+ W* ( • 04 1663974 W*  ( . 000  030544  W* (-. 00262 5734W* ( . 0005 
141  25  4Wk<  .00  0293  33-***  .00  01  3558)  )))  ) 

P S JO=FO*COS<  Z-PO) /SORT! Z> 

2 CONTINUE 
PSY0=-1 . 0F7  5 

I F ( x.FO. 0* ) GO  TO  3 
I F ( Z . GT . 3 . ) GO  TO  4 

P SYO =0.63661 977* ALOG ( 0 «5*Z ) *BSJ0  40 . 3 6 74660 1 4Y* ( 0 . 60550 3664 Y* ( - 0 

1 ,743503844Y*(  0.253  0 01  1 74Y*< -0.  042  61  21  44 Y* ( 0. 0042 7 91  6- Y*  0 . 0 00? 484 6 ) 

2 ) ) ) ) 

GO  TO  3 

4 CONTINUE 

fsyo=fo*sin(z-po)/sopt{ z) 

3 CONTINUE 
HANKO?=BSJO-U*BSYO 
RETURN 

ENO 


These  function  subprograms  compute  the  necessary  Hankel  functions 


hi2)(x)  “ Jl(x)  “ 3 Yl(x) 

and 

hJ2)(x)  - JQ(x)  - j Yq(x) 


by  using  polynomial  approximations  as  given  in  [10,  Sec.  9. A]. 
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B-4.  Description  and  Listings  of  Plot  Subroutines 

Following  is  a list  of  subroutines  used  to  plot  Figures  6 
through  21  in  this  report.  The  driver  program  is  included  here  as 
an  example  of  how  the  subroutines  are  called.  Punched  output  was 
generated  by  the  main  program  described  in  Appendix  B-2  for  use 
by  the  driver  program.  In  each  of  the  subroutines  the  initialization 
statement 

CALL  PLOTS ( 0,  0,  0,  20.,  11.) 

is  made.  This  is  standard  procedure  at  Syracuse  University  and  simply 
tells  the  operator  that  a 20"  by  11"  space  is  to  be  reserved  on 
the  plotting  roll. 


LISTING  OF  DRIVER  PROGRAM 
SJOQ  AUCK  , T I VE=2,PAGES=40 

C PROGRAM  TO  PLOT  FIGURES  6 THROUGH  21... 

DIME  NS  I ON  C UR { 30  0 ) .PHASE ( 33  0 > »FP (300)  . ANG( 300 ) , ANG  2 ( 300 ) 
DIMENSION  GA( 300) , SF ( 3 ) 

DATA  SF/1..1...5/ 

100  FORM  AT (15) 

101  FORMAT (39X.2E 17.7  ) 

102  FGRMAKSX.F  10.2) 

103  FORM  AT ( F20 . 1 .2E20 .7 ) 

104  FORM  AT (2E1  1 .4 . 1 IX  ,E1  1 . 4 ) 

C PLOT  FIG.  6 

READ(1 .100)  N 
N2=2*N 
DO  1 1=1.3 

DO  1 J = 1 .N2 

READ  ( 1 . 1 01  ) CURC  J + ( 1-1  )*N2)  , PHASE  { J-K  I - 1 ) * N2  ) 

1 CONTINUE 

CALL  PCUP(CUR,N,N2.1 .5.3) 

CALL  PPhAS(PHASE,N,N2. 3) 

C PLOT  FIG.  7 

DO  2 1=1.3 
DO  2 J = 1 . N 2 

READ ( 1 , 1 01  ) CUR( J + ( I -1  ) * N2 ) . PH  ASECJ-M  I— 1 )*N2  ) 

2 CONT1NUF 

CALL  PCURICUR. N.N2. 1 .5.3 ) 

CALL  P0HAS(PHASE.N.N2*3) 

C PLOT  FIG.  8 

READ  ( 1.100)  NFP 
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n -o 


3 


C 


4 


C 


5 


6 

C 


7 


8 

C 


10 


NFPJ=3*NFP 
DO  3 1= I .NFP3 

Ft  AD (lal 03 ) ANG(  I ) • GA(  I ) »FP(  I) 
ANG2  ( I ) = AN  G(  I ) 

CONT I NUE 

CALL  PFP(ANG,FP,NFP, 3) 

CALL  PGA ( ANG? . GA . NFP • 1 . • 3) 

PLOT  FIG.  9 

Ft  ADI  1 . 100)  NFP 
NFP3=3*NFP 
DO  4 1=1  • NF  P 3 

REAO(l.lOJ)  ANGI  I )iGA(  I)  .FP(  I) 
ANG2  I I ) = ANG I I > 

CONT  INUE. 

CALL  PFP( ANG.FP.NFP, 3) 

CALL  PGA(ANG2,GA,NFP,1 ..3> 

PL-OT  FIG.  10 

READ (1,1 00)  N 
N2  = 2 *N 
NP1  = N+  1 
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DU  5 J= 1 ,N 

READ  (1,101)  CUR ( J ) 

CONT  INUE 

DO  6 J=NP1,N2 

READ (1  .102)  CUR ( J ) 

CONT  INUE 

CALL  PCUR(CUR.N.N2. 1 .. 1 ) 

PLOT  FIG.  11 

READ (1,100)  N 
N2  = 2 AN 


NP 1 = N* 1 


DO  7 J=  1 ,N 

READ ( 1.101)  C UR ( J) 

CONT  INUE 

DO  8 J=NP1 ,N2 

READ (1,102)  CUR( J ) 

CONT  INUE 

CALL  PCURICUR.N.N2, 1 ., 1 ) 

PLOT  FIGS.  12.13.  AND  1A 

READ ( 1 , 1 00  ) N 
N2  = 2 *N 
DO  9 1=1.3 

DC  10  J= 1 , N2 

RE AD (1.1 01)  CUR ( J ) , P HA  SE ( J ) 
CONT INUE 

CALL  PPHAS( PHASE. N.N2. 1 ) 
CALL  PCUR( CUR ,N.N2. SF(  I ) . 1 ) 
CONT INUE 

PLOT  FIG.  15 

READ.N. NFP, PHI  0, OPH I 
NFP3=NFP*3 
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DO  11  1 = 1 • N F P 3 

READ (1*103)  ANG<  1 ), GA(  I ) , FP ( I ) 

ANG2 ( I ) =ANo<  I 1 
1 1 CONT  IN  U t 

CALL  PFP(ANG,FP,NFP.31 
CALL  3GA( ANG2iGAi NF P » 1 . » 3 ) 

C PLOT  FIG.  16 

READ ,N , NFP , PH  I 0 .DPMI 
NFP3=NFP  *J 
DO  12  I=l,NFP3 

READ ( 1 • 1 03  ) ANGlI  ),GA1I)  .FP(  I) 

ANG 2 1 I )= ANGl  I ) 

12  CQNTINUt 

CALL  PFP(ANG,FP,NFP, 31 
CALL  PGA ( ANG2  .GA . NFP .2  . . 3 ) 

C PLOT  FIGS.  17,16,  AND  19 

READ ( 1 . 1 00  ) N 
N2  = 2 *N 
DO  13  1=1,3 
DO  1 A J= 1 , N2 

READ (1,101)  CUP ( J ) .PHASE ( J) 

14  Cl  NT  INUE 

CALL  PPHAS ( PHASE , N, N2.  1 ) 

CALL  PCURl CUP  ,N,N2, 1 . , 1 ) 

13  CONT INOL 

C PLOT  FIG.  20 

NTP-41 
DO  lb  1=1.3 
DO  15  J= 1 , N TP 
K=JF(  I-l  ) *N  TP 

READ (1.1 04  ) ANGl K) ,A NG 2 ( K.  ) • G A ( K ) 

15  CONTINUE 

CALL  PT FANS ( ANG. ANG2 . GA. NTP. 3 t -l  ) 

C PLOT  FIG.  21 

NP=2 

DC  16  1=1, NP 

DO  16  J= l.NTP 
K = J+ ( I-l  ) * NT  P 

RE  AD (1  . 10  4)  ANGIK 1 ,A  NG  2 ( K ) ,GA(K) 

16  CONTINUE 

CALL  PT  PANS ( ANG. ANG2 .GA.NTP.NP,  1 ) 

STOP 

END 
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C LIFTING  Of  SUUROU  TINE  PCUM 

SUUKCUTINC  PCUR ( Y ,N. N2 .SC.NP > 

C PROGRAM  TO  PLOT  CURRENTS 

DIMENSION  X X( 5 1 .Y Y( 5 ) . X( 300 ) 

D I ML  NS  I ON  Y 1 ( l 50  ) . Y2  ( 1 50  ) , Y ( 300  ) 

DATA  XX/2.,7.,7.,2.,2./, YY/1  ..1..5..5.  ,1./ 

NSN=  5 
SN=SC*  2. 

XW=XX( 2 )-XX ( 1 ) 

YW  = Y Y( J l-YY ( 2 ) 

SI  =YW/SN 

DX=XV*/N 

0X2=DX/2. 

CALL  PLOTS! 0. 0 i0 .20 ., 1 1 . 1 

CALL  L I NE< XX ( 1 > . Y Y( 1 ) , 5.  1 ,0 . 0 1 

J T = 1 

NPl=N+ 1 

DO  2 J-l  .2 

YS=Y Y( I ) ♦( J- l ) *YW 

XS=XX ( 2 ) -( J-l ) *XW 

DO  1 1=1* NP 1 

XI  = XXt  l ) M J - l >*XW*(  1-1  ) *D  X*  J T 
CALL  SYMH0L  (XI, YS, .14,13,0..-!) 

1 CONTINUE 

DO  3 1=1 ,NSN 

Y I =Y  Y ( l ) H J-l  1 *YW ♦<  1-1  ) * JT 
CALL  SYMUl'L  <XS.YI..14,13,90..-I> 

3 CENT INUE 
JT=- JT 

2 CENT  INUE 

SN2=  SN/4 . 

DO  4 1 = 1 ,NSN 

Y I =Y  Y ( 4 ) — ( I — 1 1-.0  7 

CALL  NUM0ER1  XX  ( 1 )-.(>4,YI  , . 1 ♦ , SN  , 0 . , 2 ) 

SN=SN-SN2 

4 CCNT I NUb 

CALL  SYMUUL (XX(l)-.7,YY(l)M.4,.l4 , 9HM AGN I TUDE , 9 0 . .9 ) 
CALL  SYMOUL (XXI l 1-.0J.YY1 1)-. 27. .14, 112, 0..-1) 

CALL  SYMBOL  ( XX  ( 1 ) ♦XVX/a  . , YY  ( 1 ) , . 2 , 1 3,  0 . • - 1 ) 

CALL  GYM DLL I XX ( 1 ) 4XW/2 •-  • 1 7 , Y Y ( I ) - . 2 7 • . 1 4 , • * /2 • , 0 . , J ) 
CALL  SYMBOL (XX<1)42.00,YY(1)-.6..14,8HPUSITION,0.,8) 
CALL  SYMBOL (XX(2)-.0  35  .YY<1  )-.27,.14,l02,0..-l  1 
DO  5 I = 1 ,N 

X<  I ) =XX  1 l ) ♦ 1 I-  .5  1 *l)X 

5 CONT INUL 

DO  10  K=  1 , NP 
DO  (<  I = 1 ,N 

Y 1 ( I 1 - Y 1 I ♦ ( K-  1 ) *N  2 ) * ST  ♦ Y Y<  1 ) 

«•  CUNT  INUE 

'll  t I =NP1  , N2 

».  I I N>  VI  I ♦ C K — I I *N2  1 * SF  ♦VY  ( I ) 

• ( ONT  | N< II 


*2 
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OH  0 | - 1 .N 

cAl  L t>Y  MUl.l  IX(I).YIII)i.OM)iOm-I) 
CAL  L LtY  MiUlL  ( X<  I ) . YP  ( 1)  I . or , 2 .0  • ,-l  ) 
CUNT  INUf 

CAl  l l INI  (X(l)*Vltl)»N.l»OfO> 

CAll  LINI  ( X(  I ) iVcM  1 I iN  ■ 1 tOiO) 

CUNT  IN  til 
Cl  1 UKN 


This  subrout  1 no  plots  tbo  magnitude  ot  and  M,  versus 
their  position  on  the  (aee  1'^  or  1’,  and  connects  the  points  with 
straight  lines.  The  argument  parameters  are: 

V - arrav  containing  | N ^ | and  |M , | sequent  la 1 1 y stored. 
N » number  ot  expansion  functions  on  slit  face  1'. 


and  I',. 

N2  - 2*N 

SC  * scale  factor  such  that  the  maximum  value  ol  the 
plot  ordinate  equals  St'*.’. 

NT  » number  of  plots  ot  ] M ^ | and  |M,|  tv'  be  drawn  tv'i 
one  picture. 

The  arrays  having  minimum  dimensions  are  Y(N2*N1'1,  XtNl.  YltN'',  and 
Y2(N).  All  v't  the  statements  up  tv'  I'O  loop  S generate  the  picture 
frame  with  labels.  Ittslvle  00  loop  10,  the  values  v't  | M ^ | and  |M,|  are 
taken  out  v'l  arrav  Y,  placed  in  the  dummy  arrav  Y1  and  Y.',  anvt  markevl 
with  squares  and  triangles  respectively.  This  is  done  t v't  NT  arrays 
of  |M  | ami  |Mj . 
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C LIST  IN  o OF  SUBROUTINE  PPHAS 

C 

SUBROUTINE  PPHASI Y.N.N2 ,NP) 

D I ML  NS  U'N  X X ( 5 ) • Y Y ( 6 ) • X(  SO  • Y ( 300)  . XK  ( 5 ) 

DATA  XF/.5  • .35. .23, .5. .65/ 

DATA  XX/2..7..7..2..2. /, YY/l ..1..5..5. .1./ 

XW=XX( 2 ) -XX ( l ) 

Y*  = Y Y(  J)-YY ( 2 ) 

NSN  = 5 
DX=X  W/N 

CALL  PLOTS( 0,0.0.20.  .1  1.) 

CALL  L INE ( XX ( l ) , Y Y(  1 ) . 5,  1 . 0 . 0 ) 

JT=  1 
NPl=N4 1 
DO  2 J-  l .2 
Y5=YY(  1 )4( J-l  ) *YHI 
X o— X X ( 2)  — ( J-  1 ) *XW 
DO  l I=1,NP1 

X 1 =X  X ( 1 )4(  J-l  )+XW  + ( 1-1  ) * l)X*  J T 
CALL  SYMBOUXI  .YS..I  A,  13.0..-1  > 

1 CUNT  I NUt 

DO  3 1-liNSN 

YI=YY(1  ) 4 ( J-l  ) * Y W 4 ( 1-1  > * JT 
CALL  SYMBOL ( XS.YI  ..14,  13.90.  .-1) 

3 CONTINUE 
JT  =—  JT 

2 CONTINUE 
S N= 1 80  . 

DO  4 1=1 ,NSN 

YI=YY{4  )-(  1-1  ) 

CALL  NUMOL  M(XX(1)-XF(I),YI-.06,.14,SN,0.,-1) 

SN=  SN— 9 0 • 

4 CONT  INUE 

CALL  SYMBUL  <XX(1)-.7,YY(  1 ) 4 1 .67  , .1 4 .5HPHASE. 90 . , 5) 

CAl  L SYMBOL  (XX(  1 >-.0  3.  YY  ( 1)-. 27. .14, 112,0.  . — 1 ) 

CALL  SYMBOL (XX(1  ) 4XW/2 . . YY(  l ).  .2.13.0..-1 ) 

CALL  SYMBOL  < XX  ( 1 ) 4X  W/2  .-  . 1 7 • Y Y ( 1 )-.27,.14,»W/2*  ,0..3> 

CAl  L 5 Y MUOL  <XX(2)-.0  35.YY(  1 )-.2  7..14,102.0..-l) 

CALL  SYMBOL (XX(l)42..YY(l)-.6,.l 4 , 8HPOS IT  ION ,0  . . 8 ) 

DO  10  K=1.NP 
DU  S I = 1 . N 

X<  I 1 =XX<  1)  4 ( I — . f>  ) *DX 

5 CONT  1NUL 
DO  6 I = l ,N 

Y l I ) =Y ( 1 ♦ C K — 1 ) *N2 )/9  0. 43. 

6 CONT I NUL 

DO  7 I = NP 1 , N2 
Y ( I ) = Y ( I ♦ ( K — 1 ) *N2  1/90.4  3. 

7 CCNT INUL 
DO  8 J=  I »N 

CALL  SYMBOL (X(J).Y(J).  .07.0.0..-1) 


CALI  SYMBUL (X(J),Y(J4N)..07,2.0..-1) 
8 CONT  INUl 
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CALL  L I NL(X ( l ) .Y ( 1)  , N,  I ,0,0  ) 
CALL  LINLI XI  1 > .YCNPl  > ,N,  1 ,0 .0) 
10  CONT INUE 
R1  TURN 
END 


This  subroutine  plots  the  phase  of  and  M.,  versus  their 
position  on  slit  face  or  and  connects  the  points  with  straight 
lines.  The  argument  parameters  are: 

Y ■ array  containing  the  phase  of  and  M.,  in  degrees. 

N - number  of  expansion  functions  on  slit  face  and  I\, . 
N2  - 2*N 

NP  “ number  of  phase  plots  desired  per  picture. 

The  arrays  having  minimum  dimensions  are  Y(N2*NP)  and  X(N) . All  of 
the  statements  up  to  DO  loop  10  generate  the  picture  frame  with 
labels.  DO  loop  10  then  takes  the  phases  out  of  array  Y and  plots 
the  phase  of  with  squares  and  the  phase  of  M9  with  triangles. 

This  is  done  for  NP  arrays  of  and  M?. 

C LISTING  OF  SUBROUTINE  PTP 

C 

SUBPOUT  I NL  PFP ( X » Y « NFP  * NP } 

C SUBROUTINE  TO  PLUI  NORMAL  I Zt  D IILLO  PATTI  KN  • 

D I ML  NS  ION  XX(S).YY(5).X( J00 ).Y(300) 

DATA  XX/?«  . 8 . • 6 ■ * 2 • .2./.YY/1  . • 1 • • S • » S . • 1 • / 

CALL  PLOTS! 0.0.0*20. .1 1. > 

CALL  L INL( XX  I 1 ) *YY(  1 ) *b. 1 .0. 0 ) 

XK*XX( 2 >-XX(  1 ) 

YW=YY! J)-YY!2I 

nx=xw/o . 

J T = i 

NY  = 6 
NX*  7 

DO  2 J = 1 . 2 

YS=YY( 1 ) ♦( J-l ) *YW 
XS=XX<  2)-(  J-l  » *XW 
DO  1 1*1. NX 

XI  *XX!  I ) ♦(  J-l  )*XW+(  1-1  )AI)X*  JT 
CALL  SYMDUL ( XI  .Yb. • 1 4*  13*0* *-l > 

1 CONTINUE 


77 


THIS  PAGE  IS  BEST  QUALITY  FRACTICABU 
JHOM  COtY  FURNISHED  TO  DDC  — 


DO  3 1=1, NY 

Y1=YY<  l ) ♦ C J — l ) * Y W ♦ ( 1-1 
CALL  SYMUUL  < Xb . Y l i < 1 4 ■ 1 3 . 90  . . - 1 ) 

3 CENT INUC 
JT=— JT 

2 CUNT  IN UF 

bN=  1 . 

DO  A 1=1  , N Y 

YI=YY(4)-( l-l)».3-.07 

CALL  MUMULH ( XX  ( 1 ) - .5 . Y I . . 14 , SN. 0 . . 1 ) 

SN=SN-.2 

4 CONT  INUE 

CALL  SYMUOL  (XX  < 1 ) -.7  ,Y  Y(  1 ) 4 .32.  . 14, 244NURMAL  I 7t O FIELD  PATTERN, 
*90.  . 24  ) 

GN=  — 90  . 

DO  5 1=1,3 
X 1 = X X ( 1 )-.  1 7*1  1-1  ) 

CALL  NURUL  R ( X I , YY  ( 1 1 — • 3,  .14.SN.0..-1  ) 

SN=  SN+  3 0. 

5 CONT INUE 

CALL  5 Y MUOL  <XX<l)4XW/2.-.04»YY{I)-.3..14.ll2.0..-l) 

SN=30. 

DO  ft  1=5,7 
Xl=XX(  1 )—.!♦(  I - 1 ) 

CALL  NUR  JLR (XI ,YY(l)-.3..l4,SN,0..-l) 

SN=S  N4  30 . 

6 CONTINUE 

CALL  bYMUOL  (XX(  1 ) *2 . 7, YY ( 1 ) - .6 • . 14, 5H ANGLE . 0 . , 5) 

DO  11  K = 1 » NP 
DU  10  I = 1 . NF  P 

X(l)=XX(l)4<X4/180.)*(X(l4(K.-l)*NFP)490.) 

Y ( I ) = Y(  I *-  ( K—  1 ) *NKP>*  Y.V4Y  Y ( 1 ) 

10  CONTINUE 

CALL  L I NL ( X ( 1 ) , Y ( 1 1 . NF P,  1,0,0) 

1 1 CONT  I N Ut 

RCTUPN 
END 


This  subroutine  plots  the  normalized  far  field  pattern 
measured  in  the  half  space  region  c as  a function  of  the  observation 
angle  measured  from  the  x axis.  The  argument  parameters  are: 

X ■ array  containing  values  of  the  angle  at  which  the 
far  field  measurements  are  made. 

Y = array  containing  far  field  measurements  normalized 
so  that  the  largest  measurement  equals  unity. 
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C LIFTING  Of  SUBROUTINE  PGA 

C 

SUBROUTINE  PGA( X, Y ,NF^  ,SC ,NP ) 

C SUDROUI  INE  TO  PLOT  GAIN  PATTERN. 

D I ME  NS  I ON  X ( 30  0)  , Y { 3 00 ) , XX ( 4 ) . Y Y { 4 ) 

DATA  XX/2.,7.,3.,J./,YY/4.,4.,2.  .6./ 

DATA  NX/4/ . NY/S/ . PI /3.  14  I 593/ 

CALL  PLOTS! 0.0 .0 ,20.  ,1  1 • ) 

CALL  L INE(XX(  1 ).YY(  1 ),2,  1,0.0) 

SN=4  . * SC 
DO  1 1=1, NX 

X I =X  X ( 2 ) — ( I — 1 ) 

CALL  SYMBOL (X1,YY(1),.14,13,0.,-1) 

CALL  NUMOC  R (XI-.03.YY!  1 ) - . 22  . . 1 4 , SN . 0 . . - 1 ) 

SN=  S N — SC 

1 CONTINU! 

CALI  L 1 NE  (XX(3).YY(J),2.  1,0,0) 

DO  2 1=1, NY 

Y I = Y Y ( 4 ) — ( I — 1) 

CALL  SYMBOL  (XX(1)M..YI,.14,13,90..-1) 

2 CONTINUE 

XN=  X X ( 1 I+.H6 

CALL  NUMOC R(XN,YY(4)-.07..14,2.*SC,0..-1) 

CALL  NUMBER ( XN,YY(4)-1 .07, . 14, SC, 0., -I ) 

CALL  NU  MOL R(XN.YY(3)-.07,.14.2.*SC.0..-1) 

CALL  NUMOLR ( XN  , Y Y ( 3 ) ♦ . 93  . . 1 4 , SC . 0 . , - 1 ) 

DO  4 K = l , NP 

DO  3 1 = 1 , NF  P 

YT  = Y IO(K-l)  *NrP) 

Y(  I )=YY ( 1 ) ♦ Y ( I ♦ ( K — 1 )*NFP)*S I N ( P I *X( I ♦ ( K— 1 ) *NFP )/ 18  0.) /SC 
X(  I ) =XX(  1 )4YT*C0S(PI*X( I M K— 1)4  NFP  I / 1 30 • ) /SC4 1 . 

3 CONT  1NUE 

CALL  L I NL ( X ( I ) ,Y(  l ) . NFP.  1 , 0, 0 ) 

4 CONTINUE 

CALL  S Y MOLL (XX(3)+.5.YY( 3)-. 5, • 14, 'GAIN  PA T TCKN • . 0 . . 1 2 ) 

RE  TURN 

END 
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This  subroutine  plots  a polar  pattern  of  the  gain  computed 
in  the  main  program  for  the  half  space  region  c.  The  argument  para- 
meters are: 

X = array  containing  the  gain  measurement  angle 
measured  from  the  x axis. 

Y = array  containing  the  gain  measurement. 

NFP  = number  of  gain  measurements  made  per  pattern. 

SC  = scale  factor  such  that  A _<  SC*maximum  gain  value. 
NP  = number  of  patterns  desired  per  picture. 

The  arrays  with  minimum  dimensions  are  X(NFP*NP)  and  Y(NFP*NP). 

All  of  the  statements  up  to  DO  loop  4 generate  the  axes  with  labels. 
The  data  is  then  scaled  and  plotted  using  polar  coordinates. 
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C LISTING  or  SUBROUTINE  PTRANS 

C 

SUUROUT INC  PTRANS (X.Y. T.NTP.NP. ID) 

DIMENSION  X { 300 > . Y( 3 00 ) . T< 300 ) 

DIMENSION  XXI  5)  * Y Y < 5 ) 

DIMENSION  DUMX <8o ) , OUMYI 86 1 

D I ML  NS  I ON  TMF 16) , PHI  (8 ) , TMH I (8)  .TMDI8 ) t PHD (8 ) ,TMol  (8 ) 
OAT  A PH  I/O  . . 1 0 . .20.  . 30  . . A0.  , 50.  . 60.  . 70.  / 

DATA  TMF/1  . , .975. .025.  .86, .75,  .65,  .56.  .485/ 

DATA  TMf  1/.  98,.  95,  .92.  .85,.  7«J  , . 725,  .66,  .62/ 

DATA  X X/2  « »6.,6.»2.»2./ » Y Y/ 1 . • l • ,4. ,4.,  1 ./ 

DATA  PI/3.141593/ 

XW=XX( 2 ) —XX ( I ) 

Y*=YY( 3 )-YY<2 ) 

DX=XV*/9  . 

DY=YW/1 0. 

CALL  PLE TS(  0 ,0 .0 ,20.  .1  1 . ) 

CALL  L INEIXXI  I >.YY(  1 ),5.  1,0.0) 

JT=1 

DO  2 J= 1 ,2 
YS=YY( 1 )♦( J-l ) * YW 
X S=X  X I 2 ) — ( J-l  )*XW 
DO  1 1=1,10 

XI=XX<  l ) ♦ ( J — 1 )*XIKt(  1-1  ) * DX  * J T 
CALL  SYMQOL ( XI  ,YS » . 1 4 . 1 3 .0.  .-1 ) 

1 CONT  INUE 
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DO  3 1=1,11 

YI=YY<  1 )♦(  J-l  ) * YW  4-  ( 1-1  )*  JT*DY 
CALL  S Y MBOL ( XS  • Y I . . 1 4 , 1 3 . 90 . , - 1 ) 

3 CONTINUF 
JT=- JT 

2 CONTINUE 

SN=1  • 

DO  4 1=1,3 

YI=YY( 4 )-<  1-1  )* 1 .5-.07 

CALL  NUMBE  R<  XX  ( 1 ) - . 50  , Y I • • 1 4 , SN  , 0.  . 1 ) 

S N=  SN-.b 

4 CONTINUE 

CALL  SYMBOL  ( XX < 1 > - . 6 , Y Y ( 1)4-1. 13,. 14, 'T*  ,90.  .1) 

CALL  SYMBOL  (XX(1  J-.6.YYI  1 ) 41  .33,  .14,  'CUS  • , 90  . , 3 ) 

CALL  SYMBOL  ( XX  ( 1)-.6,YY<1>41.77,.14,36,90..-1) 

CALL  NUMBER (XX ( 1 ) - . 0 35 , Y Y ( 1 » - .25  . . 1 4 , 0 . . 0 . . - 1 ) 

CALL  NUMBLRt  XX(2)-.11S,YY(1  )-.25,.14,90.,0.,-1) 

CALL  SYMBOL  ( XX(  1 ) 4-.  75,  YY  ( l )- . 4,  . 14  , • ANGLE  OF  I NC  I DC  NCE  • , 0 . , 1 8 ) 
DO  6 J= 1 , NP 
DO  5 I = 1 ,N  TP 
K=  I 4-  ( J — 1 ) *N  TP 

DUMX ( I )=XXC  2)-{ 2.4XW/P I ) * ( AOS (AT AN2 ( Y ( K),X(K)))-PI/2.) 

DUMY  ( 1 >=T(  K ) *YW4Y  Y(  2 ) 

5 CONTINUE 

CALL  L I NE ( O UMX ( 1 ) ,DUMY{ 1),NTP, 1,0,0) 

6 CONT  INUE 

IF ( I D.EQ.- 1 ) GO  TO  10 
DO  8 1=1,8 

PHD ( 1 ) =XX( 1 ) 4XW*P H I ( I ) /9  0. 

TMO<  I )=YY(  1 ) 4-  T MF  ( I)*YW 
T MO  1 ( I )=YY(  1 ) 4TMF  1 ( I ) *YK 

CALL  SYMBOL (PHD( I l.TMOll l),.07,2,0..-l ) 

CALL  S Y MBOL ( PHD ( I ).TMO(I)  ,.07,0,0..-!) 

8 CONTINUE 

10  CONTINUE 

RETURN 
END 


This  subroutine  plots  Tcosiji  versus  $ where  T is  the  transmission 
coefficient  computed  in  the  main  program  and  <{>  is  the  angle  of  incidence 
measured  from  the  negative  x axis.  The  argument  parameters  are: 

X = array  containing  x coordinate  of  line  source. 

Y = array  containing  y coordinate  of  line  source. 

T ■ array  containing  values  of  Tcos^  where  41  “ tan  ^1^1 
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NTP  ■ number  of  values  of  T computed. 


NP 

ID 


number  of  plots  of  Tcos<t>  versus  4>  desired  per  picture. 
Integer  option  variable.  If  ID  ■ 1,  the  sample  data 
placed  in  arrays  TMF,  TMF1,  and  PHI  will  be  plotted 
with  symbols, i.e. , for  comparison  purposes.  If  ID  ■ -1, 


this  data  is  ignored. 

The  arrays  with  minimum  dimensions  are  X(NTP*NP) , Y(NTP*NP),  T(NTP*NP) , 
DUMX(NTP) , and  DUHY(NTP).  All  of  the  statements  up  to  DO  loop  6 
generate  the  picture  frame  with  labels.  Each  angle  <J>  is  computed  in 
DO  loop  5 and  stored  in  the  dummy  array  DUMX.  Array  T is  scaled  and 
stored  in  array  DUMY.  Straight  lines  are  then  drawn  between  consecutive 
points. 


I 
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